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ABSTRACT 

The rapid advance of gravitational-wave (GW) detector facilities makes it very important to estimate the 
event rates of possible detection candidates. We consider an additional possibility of GW bursts produced 
during unbound orbits of stellar mass compact objects. We estimate the rate of successful detections for spe- 
cific detectors: the initial Laser Interferometric Gravitational-Wave Observatory (InLIGO), the French-Italian 
gravitational-wave antenna VIRGO, the near-future Advanced-LIGO (AdLIGO), the space-based Laser Inter- 
ferometric Space Antenna (LISA), and the Next Generation LISA (NGLISA). The dominant contribution among 
unbound orbits that have GW frequencies in the sensitive band of the detectors correspond to near-parabolic 
encounters (PEs) within globular clusters (GCs). Simple GC models are constructed to account for the compact 
object mass function, mass segregation, number density distribution, and velocity distribution. We calculate 
encounters both classically and account for general relativistic corrections by extrapolating the results for infi- 
nite mass ratios. We also include the cosmological redshift of waveforms and event rates. We find that typical 
PEs with masses m\=m.2= 40M© are detectable with matched filtering over a signal to noise ratio S/N = 5 
within a distance d h ~ 200 Mpc for InLIGO and VIRGO, z = 1 for AdLIGO, 0.4 Mpc for LISA, and 1 Gpc 
for NGLISA. We estimate single datastream detection rates of 5.5 x 10~ 5 yr _I for InLIGO, 7.2 x 10~ 5 yr _I for 
VIRGO, 0.063 yr" 1 for AdLIGO, 2.9 x lO^yr" 1 for USA, and l.Oyr" 1 for NGLISA, for reasonably conserva- 
tive assumptions. These estimates are subject to uncertainties in the GC parameters, most importantly the total 
number and mass-distribution of black holes (BHs) in the cluster core. In reasonably optimistic cases, we get 
> 1 detections for AdLIGO per year. We can expect that a coincident analysis using multiple detectors and 
accounting for GW recoil capture significantly increases the detection rates. We give ready-to-use formulas to 
recalculate the estimates when these input parameters become better-determined. In addition, we provide the 
partial detection rates for various masses. The regular detection of GWs during PEs would provide a unique 
observational probe for constraining the stellar BH mass function of dense clusters. 
Subject headings: gravitational waves - black holes - globular clusters: general 



1. INTRODUCTION 

Interferometric gravitational-wave (GW) detectors LIGO, 
GEO, TAMA, and VIRGO are searching for GW signals with 
unprecedented sensitivity (Hughes et al. 2001; LIGO Scien- 
tific Collaboration 2005a; Abbott et. al. 2005a,b; Grote et 
al. 2005; Ando et al. 2005; Acernese et al. 2005). For LIGO, 
the noise levels are already reaching the goal level necessary 
for the detection of the strongest signals. It is very impor- 
tant to analyze the detection capabilities of these detectors 
and to estimate the rates of potentially detectable GW sig- 
nals. There is already a considerable list of possible detection 
candidates (for a review see Cutler & Thorne 2002): the inspi- 
ral of neutron star (NS) or black hole (BH) binaries, the tidal 
disruption of NS by BH in NS-BH binaries, BH-BH merger 
and ringdown, low-mass X-ray binaries, pulsars, centrifugally 
hung-up proto neutron stars in white dwarf accretion-induced 
collapse, supernova core collapse, gamma ray bursts, and the 
stochastic background. In this paper we consider an addi- 
tional possibility, GWs produced by unbound orbits. As we 
will show, among unbound orbits near-parabolic encounters 
(PEs) produce gravitational radiation with typical frequencies 



appropriate for detection with terrestrial facilities. For close 
PEs the gravitational radiation is short and intensive, that is 
observable to large distances. Here, we estimate the expected 
event rate of detections for specific current and near-future 
GW detectors. 

Initial order-of-magnitude estimates on the detectability of 
GWs emitted during scattering and near collisions of stellar 
mass compact objects in active galactic nuclei and globular 
clusters (GCs) were made by Dymnikova, Popov, & Zentsova 
(1982). Although their study primarily focused on BH-star, 
star-star encounters, and did not provide numbers for BH- 
BH encounters, they identified these encounters to be "quite 
rare". However, Dymnikova, Popov, & Zentsova (1982) used 
an overly simplified GC model in which the velocities and 
masses of all objects were identical, and the spatial distribu- 
tion was assumed to be homogeneous. We extend the detec- 
tion rate estimates to account for the stellar BH mass func- 
tion, mass segregation, and mass-dependent relative veloci- 
ties. We show that this improvement significantly increases 
the event rate, by approximately a factor of 10 2 . In addition, 
interferometric GW detector technology has improved greatly 
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and detailed sensitivity curves are now available. Dymnikova, 
Popov, & Zentsova (1982) estimated a maximum visible dis- 
tance of D max = 20Mpc, which is a factor of <~ 100 less than 
Advanced LIGO's (AdLIGO) capabilities (see Fig. 4 below), 
i.e. a factor of <~ 10 6 less in the accessible volume of sources. 
Combining these factors, our detection rate estimates yield 
~ 10 8 times larger results for AdLIGO. 

Gravity waves emitted during PEs are also important for 
creating relativistic orbits by gravitational radiation reaction 
around the supermassive black holes (SMBH) in the centers 
of galaxies. The GWs emitted by the later inspiral of the star 
or compact object around the SMBH are possibly detectable 
by the future space detector LISA (Sigurdsson & Rees 1997; 
Freitag 2003; Gair et al. 2004) and also by ground-based 
detectors for highly eccentric orbits (Hopman & Alexander 
2005). In the present paper, we do not consider encounters 
with SMBHs, but focus on the direct detection of GWs from 
unbound orbits of two stellar mass compact objects (COs). 

Stellar mass unbound orbit encounters are expected to be 
most likely from dense star clusters with a large fraction of 
COs. Among regular star systems, these features are car- 
ried by galactic nuclei and GCs, where central densities reach 
10 4 - 10 7 pc" 3 within a region of 0.5 - 3pc (Pryor & Mey- 
lan 1993), the inner regions contain a CO fraction of q > 1/2 
(Sigurdsson & Phinney 1995). In the present paper, we esti- 
mate PE event rates for GCs. 

As compared to other GW burst sources, the big advan- 
tage in detecting PE events is that the possible signal wave- 
forms are much more reliable as the physics behind them 
is well understood. The waveforms are known analytically 
for the case of arbitrary masses moving with arbitrary veloc- 
ities but at small deflection angles (often referred to as grav- 
itational bremsstrahlung, see Kovacs & Thorne 1978), arbi- 
trary unbound orbits but low velocities in the Newtonian ap- 
proximation (Turner 1977), in the post-Newtonian approxi- 
mation (PN, Blanchet & Schafer 1989, including corrections 
C(v 2 )), in the 2PN approximation (C(v 4 ) Blanchet et al. 1995; 
Mikoczi, Vasuth, & Gergely 2005), and most recently, in the 
3PN approximation (Blanchet et al. 2005, C(v 6 )), and the 
exact numerical solution is available for extreme mass ra- 
tios using a Schwarzschild background approximation (Mar- 
tel 2004), and finally, for head-on collisions with large veloc- 
ities (D'Eath & Payne 1992). Thus, PE waveforms are avail- 
able for a very large portion of the parameter space. Wave- 
form templates can be constructed a priori, similar to inspi- 
rals. The prior knowledge of the possible waveforms allows 
the method of matched filtering detection, which helps to re- 
duce the minimum signal-to-noise ratio necessary for a con- 
firmed detection (Flanagan & Hughes 1998). 

This paper is organized as follows. In § 2, we summarize 
the relevant characteristics of interferometric GW detectors. 
In § 3, we review the PE waveforms that we adopt. In § 4, we 
describe the population models that are necessary to estimate 
the scalings of parameters and the overall PE event rates. In 
§ 5, we derive the expected number of PE event rates, cal- 
culate their maximum distance of detection, and estimate the 
implied rates of successful detections. Finally, in § 6 we sum- 
marize our conclusions and in § 7 discuss the limitations and 
the implications of this work. 

2. OVERVIEW OF GRAVITATIONAL- WAVE DETECTORS 

The new generation of GW detectors rely on interferometric 
monitoring of the relative (differential) separation of mirrors, 
which play the role of test masses, responding to space-time 
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FIG. 1 . — Goal sensitivity curves for interferometric GW detector facilities: 
InLIGO, VIRGO, AdLIGO, LISA, and NGLISA. 



distortions induced by the GWs as they traverse the detectors. 
The effect of a GW is to produce a strain in space, which 
displaces the mirrors at the ends of the arms by an amount 
proportional to the effective arm length and GW strain. For 
GWs incident normal to the plane of the detector, and polar- 
ized along the arms of the detector, the mirrors at the ends of 
the two arms experience differential motion. Waves incident 
from other directions and/or polarizations also induce differ- 
ential motion, albeit at a smaller level. 

Presently, there is an operational international network 
of first generation interferometric GW detectors: InLIGO, 
VIRGO, TAMA, and GEO (see § 1 for references). The de- 
sign of advanced terrestrial GW detector AdLIGO and space 
detector LISA is well on the way. There are also plans for a 
new generation of low-frequency underground detectors espe- 
cially sensitive for lower frequencies (DeSalvo 2004), which 
might be especially sensitive to PEs, which we will discuss in 
a separate paper in detail. Finally, there are plans for possi- 
ble future improvements of space detectors: Decihertz Inter- 
ferometric Gravitational-wave Observatory (DECIGO) (Seto, 
Kawamura, & Nakamura 2001), Advanced Laser Interfer- 
ometer Antenna (ALIA) and the Big Bang Observer (BBO) 
(Crowder & Cornish 2005). Their sensitivities, detection fre- 
quency bands and capabilities are quite different. For our pur- 
poses, a good approximation is to use: (1.) the InLIGO and 
VIRGO sensitivity goal (nearly reached) to assess present ca- 
pabilities; (2.) the AdLIGO sensitivity goal to assess future 
capabilities of ground based detectors; (3.) the LISA sensitiv- 
ity goal to assess future capabilities of initial space based de- 
tectors; and (4.) the Next Generation LISA sensitivity goal to 
assess the capabilities of possible further extensions to space 
detectors. 

The goal RMS noise density per frequency interval for the 
various detectors, including instrumental and confusion noise, 
is plotted on Figure 1. For LIGO we adopt Abbott et. al. 
(2004); Lazzarini & Weiss (1996), for VIRGO we adopt Ac- 
ernese et al. (2005) but for simplicity discard the narrow fea- 
tures 1 , for AdLIGO we adopt the noise estimates from its 
website 2 , and for LISA, we utilize the online sensitivity curve 

1 http : //www. virgo.infn.it/senscurve/ 

2 http : //www. ligo.caltech.edu/advLIGO/scripts/ref_des.shtml 
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generator 3 for the instrumental noise, and adopt the confu- 
sion noise estimate from Barack & Cutler (2004a). The noise 
levels of possible extensions to LISA named the "Next Gen- 
eration LISA" (NGLISA) are also provided by the sensitivity 
curve generator, which we also include in all of our calcula- 
tions. This is very similar to the planned sensitivity curve of 
DECIGO (Seto, Kawamura, & Nakamura 2001, see the more 
conservative case therein) and is just halfway between ALIA 
and BBO (Crowder & Cornish 2005, a factor of 3 difference 
in sensitivity from both). 

3. PARABOLIC ENCOUNTER WAVEFORMS 

The GW signal waveform for PE is available in a wide 
range of approximations (see § 1). We adopt Turner (1977) 
for the angular averaged waveforms, for which the interact- 
ing masses travel on classical Newtonian trajectories and emit 
quadrupole radiation. Other features such as spin-spin, spin- 
orbit interactions, and gravitational recoil, etc., are higher or- 
der perturbations which carry only a small total signal power 
in typical cases. Therefore for the sake of calculating the 
signal-to-noise ratio, it is a sufficient first-order approxima- 
tion to use these waveforms. 

Illustrative examples of PE waveforms can be found in 
Turner (1977) (see Fig. 4 and 7 therein). The waveforms 
are generally constituted of a large amplitude single peak or 
a jump in the time domain with characteristic time scale to, 
related to the relative angular velocity at the minimum sepa- 
ration ljq = vo/bo = 1 /to. Here vo is the relative velocity at the 
closest point, and bo is the corresponding minimum separa- 
tion. Turner (1977) provides a closed analytical formula for 
the total GW radiation energy spectrum dE/df. The spectrum 
is wide-band, for parabolic orbits it is zero at / = 0, it has a 
maximum near fo = ujo/2tt and a half-width ~ 1.5 fo. 

The characteristic signal amplitude is obtained from the 
GW energy spectrum as (Thorne 1987; Flanagan & Hughes 
1998) 



V3G 1 / 2 1+z 1 /dE 
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where z is the redshift, d^{z) is the luminosity distance (Eisen- 
stein 1997), and dE/df [(1 +z)f] is the total GW emitted en- 
ergy of the source at the emitted frequency. The orientation 
averaged signal-to-noise ratio is 



(2) 



where h(f) is the characteristic signal amplitude (1) and S n (f) 
is the one-sided spectral noise density (see § 2 for references 
for the particular detectors). Note that the relation (2) refers 
to an angle-averaged SNR obtained from the cube-root of an 
average of cubed signal amplitudes over different possible ori- 
entations of the source and interferometer. Since event rates 
roughly scale with volume, i.e. distance cubed or (S/N)~ 3 , 
this prescription is useful for estimating event rates (Thorne 
1987). For signals with optimal orientations, the coefficient 
4/5 in (2) is changed to 4. Note furthermore, that the 4/5 fac- 
tor is applicable for the detection rate using a single interfero- 
metric GW detector. There are already 4 interferometric GW 
detectors on Earth (see § 2), and it is possible that there will be 
a lot more in the future. A coincident analysis with multiple 
detectors can be used to improve the efficiency by increasing 

3 http : //www. srl.caltech.edu/ ~ shane/ sensitivity/ 



the total signal-to-noise ratio and also by insuring that at least 
one detector is close to the optimal orientation (Jaronowski et 
al. 1996). For this reason, the coefficient 4/5 in (2) is most 
likely pessimistic. For 1 detector in the optimal orientation 
and K — 1 identical detectors in random orientations a quick 
scaling of the coefficient is ^ 4+(4/5)(/T- 1) 1 / 2 . On the other 
hand, a relatively large SNR might be required to keep the 
false alarm rate at a sufficiently low level. For a conservative 
estimate on the PE rate we do not modify the 4/5 factor in the 
definition of SNR (2) and evaluate results for SNR = 5. 

The PE waveforms can be obtained from eq. (1) by sub- 
stituting the E(f) relationship specific for PEs using Turner 
(1977): 



h(f,fo) = 
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VF(f/foz) (3) 



2V3 c 4 dh f 

where M z = (1 + z){ni\m2) i ^ / (ni\ + ni2) l l 5 is the redshifted 
chirp mass if m\ and m.2 are the masses of the interacting ob- 
jects, fo z = /o/(l +z) is the redshifted characteristic frequency, 
f denotes the observed GW frequency, and F(x) is the Turner 
(1977) normalized dimensionless energy spectrum for dimen- 
sionless frequency x = f / fo z , for which J Q F(x)dx= 1. 

Equation (3) is the leading order (i.e. Newtonian) approx- 
imation to the GW waveform, h(f,fo). A remarkably advan- 
tageous feature of the waveform in this approximation is that 
it depends on only a single combination of the orbital param- 
eters, fo z . Although we need only utilize this form (3), we 
briefly note that it is also possible to express the waveform 
with the separation at closest point bo in the center of mass 
frame (Turner 1977): 



hif,b )- 
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tVW/o^o)]. (4) 
d L bo J 



Here fo z (bo) is the redshifted characteristic frequency for fixed 
initial velocity as a function of bo. The GW amplitude spec- 
trum h(f, bo) is roughly flat at low frequencies / < fo z , and 
decreases for higher frequencies. Equation (4) shows that 
h(f, bo) scales with b^ 1 for frequencies larger than the cutoff 

Modifications are necessary for relativistic encounters. The 
relativistic gravitational radiation waveforms and energy out- 
put has been calculated by Martel (2004) in the quadrupole 
approximation for a test particle approaching a Schwarzschild 
black hole from infinity on a quasi-parabolic geodesic 4 . In 
case the periastron distance bo is close to the unstable circular 
orbit, the GW energy is significantly increased (it has a log- 
arithmic singularity at A = 2, where A = bo/RsH and Rs H is 
the total Schwarzschild radius, and a factor of ~ 10 increase 
for A = 2.01 or a factor of 2 for A = 3). We adopt the fitting 
formula of Gair, Kennefick, & Larson (2005), which is cor- 
rect within 0.1% for orbits that avoid a collision, and scale the 
amplitudes of the Turner (1977) waveforms (3) according to 
the increase of the GW energy, E le i(X)/E m (\). 

KfJo) - -^!^—— — f v/fe 

Note once again, that the dimensionless periastron distance, 
A, is uniquely specified by the characteristic frequency, fo. We 

4 We continue to denote as "quasi-parabolic" or simply "parabolic" en- 
counters that have asymptotically zero velocity at infinity. Note that the tra- 
jectories are generally quite different from parabolas in the highly relativistic 
regime (for illustration, see Martel 2004). 
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derive explicit formulae for fo(X) in the non-relativistic New- 
tonian and relativistic geodesic approximations, when consid- 
ering the dynamics of PEs in § 5.1 and § 5.2 below. Quite 
remarkably, the characteristic frequency of the waveform is 
unchanged for relativistic orbits. Therefore, we do not ex- 
plore the effects of relativistic modifications in the shape of 
the GW signal waveform, we restrict only to correcting the 
amplitude. 

4. POPULATION MODELS 

The major contribution to the PE event rate is expected from 
dense star clusters with a large CO fraction. The most impor- 
tant systems carrying these properties are possibly GCs and 
galactic nuclei. In this paper we focus on GCs, but simple 
analytical scaling of the results allows a straightforward ex- 
tension to galactic nuclei or any other population of spherical 
star systems. 

The spatial distribution of GCs exactly traces the distribu- 
tion of galactic halos (Chandar, Whitmore, & Lee 2004) in 
the local universe. In this section we summarize the galaxy 
distribution and the GC abundance per galaxy, and describe 
the GC models which we adopt for the paper. 

4.1. Galaxy Distribution 

We utilize the local distribution of galaxies (Tully 1988). 
The accumulated number of galaxies to distance D can be well 
approximated by 

{N\ (D I Mpc) 09 for D < 3 Mpc 
N 2 CD/3 Mpc) 1 5 for 3Mpc < D < 16Mpc 
N 3 (D/ 16 Mpc) 2 4 for 16 Mpc < D < 60 Mpc ' 
N 4 (D/60 Mpc) 3 for D > 60 Mpc 

(6) 

where Ni, N 2 , N 3 , and N A are 23, 62, 1 100, and 26000, respec- 
tively. The average density of faraway galaxies is 0.03 Mpc -3 , 
but the local galaxy abundance is much denser than average. 
In Eq. (6) N ga \D) has a 45% jump at the Virgo cluster at 
D= 16Mpc. 

4.2. Globular Cluster Abundance 

Following Portegies Zwart & McMillan (2000), we adopt 
n gc = 2.9 Mpc -3 for the average GC abundance in the uni- 
verse. We roughly account for the clustering of GCs in the 
local universe by assuming that the distribution of GCs follow 
the abundance of galaxies. This assumption is consistent with 
observations (Goudfrooij et al. 2003; Chandar, Whitmore, & 
Lee 2004) suggesting that the population of GCs represent a 
universal, old halo population that is present around all galax- 
ies. The number of GCs within a distance D is then 

N gc (D) = y gc N g ' dl (D), (7) 
where y gc is a scaling constant relating the abundance of GCs 
and galaxies. Using the large scale average n gc = 2.9 Mpc -3 
(Portegies Zwart & McMillan 2000) and « gal = 0.029 Mpc" 3 
eq. (6) for D > 60 Mpc we get y gc = 100. Alternatively, y gc 
can be interpreted as the number of GCs per galaxy aver- 
aged over all morphological types. Concerning PE detection 
rates, we shall show in § 5.3 that typical observation dis- 
tances for terrestrial detectors are larger than the clustering 
scale (6). Therefore, the results are sensitive to mainly the av- 
erage abundance and are only slightly increased by the local 
clustering of GCs. 

The value n gc = 2.9 Mpc -3 for the average abundance is a 
conservative assumption. In their quick estimate Portegies 



Zwart & McMillan (2000) derived this value by adding up the 
contribution of galaxies of morphological types Sab, E-S0, 
and blue ellipticals. Recently, 12 nearby edge on spiral galax- 
ies were examined, resulting in much larger numbers, reach- 
ing <~ 1000- 1300 GCs for these particular galaxies (Goud- 
frooij et al. 2003; Chandar, Whitmore, & Lee 2004). In addi- 
tion to the morphological types considered in Portegies Zwart 
& McMillan (2000), dwarf elliptical (dE) galaxies also con- 
tribute to the overall GC numbers (van den Bergh 2005). The 
GC content of 69 dwarf elliptical (dE) galaxies have been esti- 
mated to host about a dozen GCs per dE galaxy (Lotz, Miller, 
& Fergusson 2004). Therefore, our results on detection rates 
correspond only to lower limits, which has to be scaled lin- 
early with n gc when more detailed estimates become avail- 
able. 

The GC distribution given by (6-7) is only valid for sub- 
cosmological scales. Assuming that D denotes the luminosity 
distance, <f L , in eqs. (6-7) which is a direct observable using 
the GW amplitude, the change in the cosmological volume 
element decreases the average density. 

N gc (d L ) = 4tt / — — — -y gc n sal dz (8) 
Jo ozdil 

We adopt the cosmological volume element (Eisenstein 1997) 
for a ACDM cosmology with {Q, m , ft A , h) = (0.3, 0.7, 0.7) con- 
sistent with recent observations of the Wilkinson Microwave 
Anisotropy Probe and the Sloan Digital Sky Survey (Tegmark 
et al. 2004). We find that the uncorrected volume element 
dl&Sl&dL is reduced by a factor of 0.7 for z = 0.1 (d^ = 
0.5 Gpc) and by a factor of 0.053 forz = 1 (d L = 7Gpc). Since 
GCs are believed to represent an old halo population in galax- 
ies, we do not account for additional possible cosmological 
evolution of the comoving GC abundance. 

4.3. Globular Cluster Models 

Here we define the GC models that we use to obtain typical 
PE parameters and event rates. First we define the common 
features that are the same for both of our models. We as- 
sume a total of N tal stellar mass stars, m star = M , spherically 
distributed within a typical radius R gc . Also within the clus- 
ter, is a CO population consisting of Nco = qN lot objects with 
q <C 1 that move in the background gravitational potential of 
the stars. We define the "typical CO mass" as (mco) = 10M Q . 
We adopt typical values of N iot = 10 6 , R gc = lpc, N C0 = 10 3 , 
with q = 10" 3 , and (m) ~ M (Djorgovski & Meylan 1994; 
Portegies Zwart & McMillan 2000; Miller 2002). 

We construct two different models for the distribution of 
mass, spatial coordinate, and velocities of stars. Model I is 
a simple plausible model to get the scaling of PE event rates 
on different cluster parameters. Here we assume a homoge- 
neous spherical distribution, and the COs within the cluster 
have the same mass and magnitudes of velocities. In Model 
II, we refine the assumptions to account for the distributions of 
masses, mass segregation, and relative velocity distributions. 

By comparing Models I and II, we find that while Model I 
gives the correct R gc , N tot , Nqo, and fo scalings, it underesti- 
mates the total event rate of a single GC by 2 orders of mag- 
nitudes. The comparison of Models I and II is necessary (i) to 
understand the origin of the large increase in PE event rates as 
compared to Dymnikova, Popov, & Zentsova (1982), and (ii) 
to understand the impact on the BH mass function of GCs on 
PE event rates. It is therefore emphasized that the GC model 
assumptions have a crucial importance when determining the 
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PE event rate. We believe that our Model II includes all of the 
essential features of GCs for the assessment of PE rates, and 
uncertainties are comprised of the uncertainties in the model 
parameters rather than additional fundamental processes 5 . 

4.3.1. Model I 

For the most simple model, the spatial distribution is as- 
sumed to be uniform within a characteristic radius R gc , and 
it is assumed that all COs have the same mass m C o = 10M Q 
and magnitude of velocity v v ; r . The orientation of velocities is 
isotropic, implying a velocity dispersion of a ~ v v ir- For this 
simple model, we estimate the relative velocity distribution 
of COs to be the same as the individual velocities v re i ~ v v ; r . 
The characteristic velocities can be obtained from the virial 
theorem for a uniform distribution of mass 6 : 



l3GN tot (m) 



R. 



(9) 



The PE event rate is calculated as the rate of scattering on a 
fixed target lattice, with incident velocity v v ; r . 

4.3.2. Model II 

We improve Model I with the following factors. The va- 
lidity and motivation of these assumptions is discussed below 
the list. 

1. We assume an equal number of NSs and BHs N^s = 
N B h- We introduce g B ii(m) = 0.5[ln(m max /m m i n )] _1 m" 1 
for the fractional distribution of BH masses among 
COs with mass m per mass interval dm, in the range 
m m i n < m < m max , with J g Ba (m)dm = 0.5. We de- 
fine the NS mass distribution g^s(m) as a Gaussian 
distribution with norm 0.5, mean 1.35M Q , and vari- 
ance O.IMq. Finally we define the CO distribution 
by gcoim) = gBH(m) + gm(m), which has a unit norm. 
For definiteness, we take A^h = 500, m m i„ = 5M Q , 
and m max = 6OM , implying that (otbh) — 2OM and 
(wco) — 10M©, but also calculate detection rates for 
more general BH mass distributions. The PE detec- 
tion rates are practically independent of the actual total 
number and distribution of NSs. 

2. We account for the mass segregation by assuming ther- 
mal equipartition among COs, objects with mass m 
have a velocity v m = (m/ (m))~ 1,/2 v v j r , and are confined 
within a radius, R m = (m/ {m))~ l l 2 R gc , while regular 
stars are distributed uniformly within a sphere of radius 
R gc , as in Model I. Since N 3> Nco, the background 
gravitational potential is determined by regular stars. 
See Fig. 2 for an illustration. For core collapsed models 
a modified scaling is necessary (see text below). 

3. The relative velocity for COs with masses m\ and m.2 is 



assumed to be v re i = w = [(w/ +m 2 1 )(m}] 



1/2,, 



The BH mass distribution is crucial for the analysis, since sig- 
nal rates scale with m 19 / 3 (see § 5 and the appendix below). 

5 For example, details like bound binary populations and interactions are 
not essential for PEs, and Spitzer instability and core collapse can be ac- 
counted for by choosing our model parameters accordingly (see § 7.3.3 for a 
discussion). 

6 There is a similar result for a spherical star system of polytropic distribu- 
tion with a root mean square radius R. The only difference is in the 3/5 factor 
of Eq. 9, which becomes 1/2 in that case (e.g. Saslaw 1985). 




FIG. 2. — Encounters for Model II. More massive COs mi > mi are dis- 
tributed uniformly within a sphere of smaller radius R[ < R2 and have smaller 
velocities. PEs between COs with mi and mi can take place within R\ . The 
relative velocity before the interaction is Vac = v\i(m\^mi). Note that Model 
II assumes a continuous mass function, for which N\ and N2 are in fact in- 
finitesimal. 



Unfortunately, the analysis of BH mass functions have not 
yet converged. Recent X-ray observations display evidence 
for 20 galactic BHs with masses between 4 < w/M Q < 14 
(Casares 2005), and <~ 45 ultra-luminous X-ray sources are 
identified with intermediate-mass black holes (IMBHs) with 
masses m = 10 2 - 1O 4 M (Ptak & Colbert 2004; Miller & Col- 
bert 2004; Gebhardt, Rich, & Ho 2005; Blecha et al. 2005). 
Theoretical predictions from two-dimensional simulations of 
stellar core collapse (Fryer & Kalogera 2001) lead to masses 
smaller than 20M Q with very different distributions depend- 
ing on the assumptions (fraction of explosion energy used to 
unbind the star, stellar winds, mass transfer after helium igni- 
tion). Sophisticated simulations of the initial phase of rapid 
star evolution assuming a lower metallicity for the progeni- 
tor stars (weaker stellar winds) appropriate for GCs and in- 
cluding a large fraction of binaries, collisions, and accretion 
leading to the mass buildup of BHs imply an initial smooth 
decreasing distribution of stellar-mass BHs with masses up to 
~ 60 - IOOMq (Belczynski et al. 2005) depending on model 
assumptions and cluster environments. Results are valid for 
timescales short compared to later dynamical evolution of 
the cluster. However, the final fate of the cluster remains 
highly uncertain. In small GCs, dynamical interactions of bi- 
naries might eject a significant portion of the stellar-mass BH 
population (Sigurdsson & Hernquist 1993; Portegies Zwart 
& McMillan 2000; O'Leary et al. 2006). Following Miller 
(2002) and Will (2004) we adopt gBn(m) oc m" 1 leaving the 
minimum and maximum masses free parameters. However, 
most recent population synthesis simulations (Belczynski et 
al. 2005) typically yield steeper BH mass functions. For this 
reason we felt it important to compute results for other dis- 
tributions gBH(m) oc m~ p with p = 0, 1, and 2, as well. Con- 
cerning the other parameters, in our standard model we take 
m m i n = 5M and m max = 60M Q , for which (twbh) — 20M Q and 
(mco) — 10M Q , and assume jV B h = 500, but also calculate re- 
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suits for other A^h, mmm, and m max values. 

The NS mass distribution that we utilize is supported by 
observations of 26 radio pulsars and 4 X-ray binaries (Post- 
nov & Prokhorov 2003, and references therein). The distri- 
bution is sharply peaked around 1 .35M Q . Our assumption on 
the total number of NSs (N^s = 500) is somewhat arbitrary, 
values could be higher or lower depending on what fraction 
of NSs are ejected by kicks during their formation. Our de- 
tection rate estimates can be scaled to the appropriate value 
using ^ NS -ns oc N^ s and ^ BH -ns oc N ns Nbh (see § 5 below). 
Note however, that in § 5 we show that NS interactions con- 
tribute a negligible fraction of the PE event rate (see Figs. 6 
and 7 below). 

Note that we do not consider encounters between white 
dwarfs (WDs). The detection of PEs between WDs are even 
worse than for NSs, since WD masses are even smaller (Bin- 
ney & Tremaine 1987). Our final results show that the PE 
detection rates roughly scale with ~ m 8 33 , implying that rates 
are enormously suppressed for WDs: by a factor of <~ 10 15 
relative to 50M Q BHs, which are typical for detection (see 
Fig. 6). Moreover, WDs are disrupted by tidal torques for 
close encounters, at sub-Hz frequencies. Therefore, PEs of 
WDs are completely invisible for terrestrial detectors. 

The second improvement is to account for mass segrega- 
tion. The differentiation of the stellar population with mass 
within the cluster core is a consequence of thermal equilib- 
rium (Binney & Tremaine 1987). Objects with masses larger 
than 10M Q had enough time to relax to thermal equilibrium 
within the lifetime of the GCs (Farouki & Salpeter 1982). 
This means that the kinetic energies of each of the component 
stars are drawn from the same distribution, implying that the 
typical speed of an object of mass m is v m = (mj '(m))~ 1/2 v v i r , 
causing the object to sink to the core of the cluster. For 
a nearly homogeneous distribution of background stars, this 
implies that the maximum radius available to a given mass 
is R m = (m/ (m))~ l / 2 R gc (Binney & Tremaine 1987). For the 
BH distribution given by gBH(m) we get that BHs with mass 
<~ 50M Q (which make the dominant contribution to PE rates, 
see Figs. 6 and 7 below) are confined to a radius <~ 0.14/? gc . 
Note that our scalings based on thermal equipartition might 
not hold in case Spitzer instability leads to core collapse, cre- 
ating a dynamically decoupled core of high mass BHs (typi- 
cally R C0K <~ 0.01-0. 10/? gc depending on the fraction of pri- 
mordial binaries, Heggie, Trenti, & Hut 2006). We can ac- 
count for core collapse by simply scaling our final results on 
detection rates appropriately with A^H^core^ore ( see § 7.3.3 
for a discussion). 

Finally, we discuss the assumption on the relative veloc- 
ity. The velocity distribution of stars in GCs is well described 
by the King-Michie (KM) model (Meylan 1987), which is 
roughly a Maxwell-Boltzmann (MB) distribution with a max- 
imum velocity cutoff. It is well-known that the relative ve- 
locity distribution for MB individual velocity distributions is 
also MB for the reduced mass = m 1 ni2/(m 1 + m2) (Binney 
& Tremaine 1987). Thus (v 12 ) RM s = (»*i/» 1/2 (v mi )RMS = 
[K'+m^Xm)] 1 /^™. 

Note that we do not utilize the exact velocity distribution, 
but associate the same fixed velocity value for every object 
with identical masses. Relaxing this approximation and ac- 
counting for MB velocity distributions leads to a change of 
only a few percent in the encounter rate results (the correction 
is (3/7r)'/ 2 for Model I in the range of GW detector frequen- 
cies, see Kocsis & Gaspar 2004 for the derivation). Thus, 



the velocity distribution can be safely approximated with the 
mean value. A simple explanation is the fact that GW detec- 
tor frequencies correspond to unbound encounters with nearly 
parabolic trajectories, for which the exact value of the initial 
velocity is negligible (see § 5.1 below). 

The PE event rate for component masses m\ and n%2 is cal- 
culated as the rate of scattering with incident velocity equal to 
the initial relative velocity v re i = v i 2(m 1 ,m2). 

5. PARABOLIC ENCOUNTER EVENT RATE 

We now derive the event rate for the successful detection 
of PE signals using the two models of GCs. This section is 
divided in five parts. In § 5.1 we derive the comoving event 
rate per comoving characteristic frequency bins for individual 
GCs for the two population models. In § 5.2 we derive the 
modifications necessary for relativistic encounters. In § 5.3 
we determine the signal-to-noise ratio using the specific de- 
tector sensitivity curves and determine the maximum observ- 
able distance of PEs. In § 5.4, we add up the contributions of 
all possible GCs within the visible distance and estimate the 
PE detection rates. Finally in § 5.5 we conclude the results of 
the analysis. 

5.1. Contributions of Individual Globular Clusters 
5.1.1. Outline 

When calculating the detectable event rates for specific GW 
detectors, it is desirable to express the encounter cross section 
for particular /o frequency bins. To achieve this we first com- 
pute the interaction cross-sections for given masses and initial 
orbital parameters. Then, the Newtonian equations of motion 
relate the initial orbital parameters to /q. Changing to the /o 
variable leads to the partial event rate for the given masses 
and characteristic encounter frequencies. Then, for Model I, 
it is very simple to add up the individual contributions of all 
objects within the cluster. For Model II, we utilize the spe- 
cific radial and relative velocity distributions, R m and v\2, and 
average over the CO mass distribution. 

5.1.2. Derivation of Event Rates 

The typical minimum distance between COs in these sys- 
tems is R gc /\/N ~ 10 11 km, a value several orders of magni- 
tudes larger than the typical minimum separation of an en- 
counter (10km < bo < 10 7 km for detectable frequencies, see 
Eq. [14] below). Therefore, a sufficient approximation is 
to consider short-time two-body interactions during encoun- 
ters 7 , and constant velocities in between events. The PE event 
rate can then be simply estimated by a scattering of particles 
with incident initial velocities Voo = v re i on a still target lattice. 

Since the velocities are assumed to be locally isotropic ev- 
erywhere in the cluster, the cross section of a particle with 
an impact parameter between boo and b^ + db^ is da = 
2-Kboodboa- We proceed to express the infinitesimal cross sec- 
tion for d/o bins. 

We derive encounter parameters with a non-relativistic 
Newtonian description. The separation, bo, and relative ve- 
locity, v , at periastron can be computed from the initial 
conditions of the interacting bodies. The initial parameters 
are the impact parameter b^ and the velocity, v^ = v re i = 

7 We shall discuss below that the PE event rate is not sensitive on whether 
the interacting participants are elements of regular bound binaries or if they 
are single objects. 
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Vi2(/Wi,/W2), of the scattered particle (see § 4.3.2 for the def- 
inition). Using the conservation of mechanical energy and 
angular momentum we get 



b. 



bo 



and 



: V 



(10) 



where 7 = (GM)/(Z?ovo 2 ) is the ratio of potential energy and 
double kinetic energy at the closest point 8 , with M = m\ + m 2 . 
Although Eq. (10) is strictly only valid in the comoving ref- 
erence frame of the center-of-mass, it is an adequate approxi- 
mation for the realistic parameters as long as Voo <C Vo. Again, 
we point out that the relevant encounters are nearly parabolic, 
so that the initial velocity distributions have a negligible im- 
pact on the result. 

Let us express the cross section with variables bo and Voo 
using Eq. (10) and make use of our simplifying assumption 
that Voo is a constant for fixed masses (see § 4.3). The result 
is 

GM \ 

■ + 1 )2irbodbo. (11) 



dcr = 



vLbo 



Here, the first term dominates the parenthesis for typical Voo 
and bo values (see Eq. 14. below). This term is responsible 
for the deflection of trajectories due to gravity. 

The characteristic GW frequency is directly related to the 
minimum separation, bo, and the relative velocity vq at bo by 



u} = 2nf = — . 

bo 



(12) 



Using Eq. (10), 



/ 2 u \ V 2 
.o = (2GM) 1 /V 2 (l + ^) ■ (13) 

For typical frequency bands of InLIGO (LISA) the correc- 
tion on the RHS is v 2 oc b Q /(2GM C0 ) < 10" 8 (10 -4 ) for InLIGO 
(LISA). In order to get the event rate per frequency bin, we 
need the inverse relationship bo(oJo,Voo)- This can be ob- 
tained recursively, as a power series in luq. To the second 
non-vanishing order, we get 

..2 



bo = (2GM) 1/3 w 2/3 



1 + 



1 



3 (2GM) 2 / 3 



LJ, 



-2/3 



(14) 



Substituting in Eq. (11) yields 



dcr = 



2tt (2GM) 4 / 3 



-IV, 



-2/3 



1 + 



8 Voc _ 2/3 \ d^p 
3(2GM) 2 7^ )^o~- { ) 
For Model I, the scattering rate for a single particle is 
"coVocdu, where «co is the number density of COs. Since 
there are a total of Nco particles, the contribution of all COs 
is 1/2 x Afco«coVooder, where the 1 /2 factor comes from the 
fact that the same ensemble of particles constitute both the 
targets and the injection. Thus Eq. (15) becomes 



where /100 = 100 Hz and 



v{ =(27r) -2/3A|o Wg^ffi = 6. 7 x 10-V-H17) 



3/2 3 

V ■ 

vir 



4GM C0 



1.7 x 10" 8 Hz. 



(18) 



8 7 < 1/2 for hyperbolic Newtonian trajectories. 



For Model II, the interacting objects with masses m\ and 
m 2 are distributed uniformly within radii Ri and R2 and have 
mass-dependent relative velocities = vn (see Fig. 2). Let 
N\ and N2 denote the number of particles with mass mi and 
OT2, respectively. Let us assume mi > m 2 , for which R\ < R2. 
In this case the interactions between masses m\ and m 2 take 
place only within a radius R\, where the density of particles 
with mass m 2 is n 2 =N2/(R\^ '/3). For a smooth distribution, 
N\ is the infinitesimal number of particles with masses be- 
tween mi and mi+dmi, i.e. N\ = Ncogco(m\)dmi (and ^de- 
fined similarly). The scattering rate for an injection of N\ par- 
ticles with Vqo velocities on a target density n 2 is is jVin 2 Vood<7. 
To get the total event rate for the cluster for loo bins we need 
to integrate over the mass distributions 



/>oo />oo 

dv n = / dmigco(wi) / dm 2 gco(w 2 ) x 
Jo Jo 

/ f \- 2 / 3 / f \" 4 /3' 

4(mi ,m 2 ) +V2(m u m 2 ) -f- 

V/100/ V/100 

where j/[ I (mi,m 2 ) and ^"(m^m^ are given by 

nv x ,„ ,-2/3 N co (2GM) 4 / 3 r . 2 ii 
v?(m u m 2 ) = (2ir) 2 ' v / 10 2 Q /3 

v^(m 1 , m 2 ) = (2?r) 4/ — ^ ' 



3 Rl 



dfo 
/o ' 
(19) 



(20) 



(21) 



where R> = max(/?i ,/? 2 ) and the mass dependence is implicit 
in the total mass M, Voo = V12, Ri, and R 2 (see the Appendix 
for explicit formulae). 

The mass integrals in (19) can be evaluated independent of 
the frequency, resulting in the same functional form as for 
Model 1(16). The constants for Model II are 



isf = l.9x lO-'V -1 
/f = 1.0x 10" 10 Hz. 



(22) 
(23) 



(see the Appendix for parametric formulae). 

For Model II, it is also interesting to get the relative en- 
counter rates forBH-BH, BH-NS, and NS-NS interactions. 
Integrating Eq. ( 1 9) over the corresponding mass intervals, we 
get 



^bh-bh = 0-996 ^ = 286 iA, 



1, BH-NS 

'-'i'ns-ns 



= 4.14 x \Q-^ = \.\9v\, 
= 7.04 x 10^f=0.02 v\. 



(24) 
(25) 
(26) 



The corresponding analytical formulas are given in the Ap- 
pendix. It is clear that BH-BH encounters dominate the event 
rates. The event rates of NS - NS encounters are more than 
four orders of magnitude lower! 

5.1.3. Discussion 

Therefore we conclude that Model II has a much larger 
event rate than Model I. By inspection of Eq. (19) and (20) 
the main factors responsible for this increase can be iden- 
tified. First, the CO density is increased by the CO con- 
finement in the core: n oc R m ~ 3 oc m 3 / 2 . Second, the typ- 
ical CO relative velocity inverse is increased: v^ oc m 1 / 2 . 
Third, the gravitational focusing is proportional to m 4//3 . Thus, 

v x oc (RlJRl ) x ((m 10 / 3 )/mco 3 ). The PE event rate is thus 



8 



Kocsis, Gaspar, & Marka 



highly inclined towards the high-mass end of the BH distribu- 
tion in the cluster. For a more precise treatment, the exact 
contributions of the component mass parameters are given in 
the Appendix, Eqs. (A4-A6). Note, that Dymnikova, Popov, 
& Zentsova (1982) obtained results equivalent to our Eq. (17) 
of Model I, but they focus on star - star encounters and use 
m = 4M© instead of mco = 10M Q . (Another difference is that 
they do not discuss /b-dependent differential rates, but derive 
the total PE rate based on a typical minimum separation. Us- 
ing m = 4M Q in Model I, the results are a factor of 10 3 lower 
than the rates for Model II.) 

Equations (16-18) and (22-23) give the resulting PE event 
rate per GC of all Newtonian trajectories between point 
masses for the two GC models considered. A significant 
shortcoming is that COs have finite radii and collide for suf- 
ficiently small minimum separations. Moreover the Newto- 
nian approximation breaks down for large velocities or strong 
gravitational fields. These effects are considered in § 5.2 be- 
low. 

-4/3 

Notice how small is the correction proportional to f in 
Eq. (16) for GW detector frequencies / ^> /J' 11 . Recall that in 
Eqs. (13) and (14) the expansion coefficient is proportional to 
Voo. Hence the leading order term is exact for = 0, thus 
it corresponds to parabolic trajectories. This proves our con- 
jecture that the unbound orbit encounter rate is dominated by 
near-PEs. 

In Eq. (16) the leading-order terms are proportional to 
1 / v v ir- The result is slightly counterintuitive if one identifies 
the star system with an ideal gas, since for ideal gases, the rate 
of collisions is directly proportional to Voc . In this perspective 
it seems reasonable to expect the encounter rate to be a grow- 
ing function of v^, for fixed frequency bins. The confusion 
arises from the fact that our GC models are using the oppo- 
site limit. For star systems the typical velocities are so small 
that the gravitational interaction dominates the motion of the 
stars. 9 Increasing the velocities decreases the gravitational fo- 
cusing, thereby decreasing the encounter likelihood. 

The expected rate of PE events for a single GC is plotted 
in Figure 3 for logarithmic frequency bins for the two GC 
models. The non-relativistic results presented in this section 
are plotted with dotted lines, which overlap with the rela- 
tivistic calculation below ~ 10 Hz. For higher frequency the 
minimum separation drops below <~ 6 Schwarzschild radii 
for the largest BHs and relativistic corrections become im- 
portant. Figure 3 displays that event rates are higher for 
lower characteristic frequencies, e.g. for model II for fo = 
0. 1 mHz (the minimum frequency for space detectors), we get 
1.9 x 10 _8 yr _1 GC -1 events, whereas for terrestrial detectors 
fo = 100 Hz it is only 1.9 x KT^yr'GC -1 . 

5.2. Relativistic Orbits 

Up to this point the PE event rates have been estimated 
for fixed characteristic frequencies but independent of the 
minimum separation bo and relative velocity vo. In addition 
to non-relativistic parabolic encounters, these events also in- 
clude head-on collisions, relativistic captures, relativistic fly- 
bys, and zoom-whirl orbits. Since we have used a Newto- 
nian analysis in the derivation, our results presented in § 5.1 

9 The ideal gas model is sufficient only for extremely small characteristic 
frequencies fo -C f\ (eqs. [18] and [23]), which is below the lower frequency 
limit of GW detectors. In this regime the stars' trajectories are only slightly 
deflected, implying that gravity, in terms of encounter likelihood, is negligi- 
ble. 




10 _ 8 10 _ 6 1Q _4 1Q -2 1Q o 1Q 2 
fo [Hz] 



FIG. 3. — The expected total rate of PEs produced in a single GC per log- 
arithmic frequency bin. GC Model I (solid) and Model II (dashed) results 
are shown including relativistic corrections for geodesies avoiding head-on 
collisions (see § 5.2). The dotted lines represent PE event rates in the non- 
relativistic approximation. The non-relativistic treatment is adequate for low 
frequencies for which the trajectories avoid collisions with minimum separa- 
tions of several Schwarzschild radii. Only a fraction of these events can be 
detected, depending on the distance of the GC. 



are valid for the non-relativistic parabolic encounters only. 
Here, we improve the classical calculation to account for 
general relativistic encounters of test particles moving along 
geodesies in the Schwarzschild space-time. This treatment is 
exact for extreme mass ratios, but as an approximation we ex- 
trapolate these formulas for general mass-ratios as well. 

To classify the orbits, we introduce a parameter A = 
bo/Rm, where bo is the distance at periastron, and Rsu = 
2GM/c 2 is the Schwarzschild radius of the total mass 10 . For 
Newtonian parabolic encounters, we get A = c 2 /vq. We distin- 
guish (i) non-relativistic parabolic encounters for A > 6, (ii) 
general relativistic flybys for 2.1 < A < 6, (iii) zoom-whirl 
orbits for 2 < A < 2.1, and (iv) head-on collisions for 2 < A. 
We restrict our calculations to A > 2, since this is the regime 
in which matched filtering can be carried out using the wave- 
forms of § 3. 

In this section, we improve § 5.1.2 to account for the rel- 
ativistic deviations in the trajectories. In practice, we repeat 
the derivation of § 5.1.2 to get the cross-section using the or- 
bital parameters of the geodesies of a test particle moving in a 
Schwarzschild space-time (Gair, Kennefick, & Larson 2005). 

For parabolic encounters the specific orbital angular mo- 
mentum is 

L=y/2GMb Q (l-A -1 )"' 72 , (27) 

where A = X(bo) defined above. The non-relativistic result is 
retained for A — > oo. Equating (27) to the angular momentum 
before the encounter L = frooVoo, solving for b^, and substi- 
tuting in dcr = 2-Kbocdboo we get 

GM 1-2A" 1 

This is to be compared to the non-relativistic analogue (11). 
The first term is the non-relativistic term for near-parabolic 
orbits and the A-dependent fraction describes the relativistic 

10 Here we restrict to BH-BH encounters which dominate event rates, see 
Eqs. (24-26). 
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correction. The latter decreases the cross section per unit bo. 
For A — > 2 the cross section becomes 0. For smaller impact 
parameters a head-on collision takes place, for which the pe- 
riastron distance and A is undefined. 

Repeating § 5.1.2, the next step is to change to the fo char- 
acteristic frequency variable. Since L=b\ d(f>/dr, where dr = 
(\-\~ l ) l l 2 dt is the infinitesimal proper time element along 
the geodesic at the closest approach (e.g. Misner, Thorne, & 
Wheeler 1973), from (27) we get 



d0 
~dt 



= uj q = It: f Q = (2GM) 1/2 b Q 3/2 . 



(29) 



Quite remarkably, this is identical to the result of the non- 
relativistic calculation for parabolic orbits (13). The GW 
waveforms have a peak at an angular frequency ojq for the 
non-relativistic encounters (Turner 1977, and see § 3 above). 
For relativistic zoom-whirl orbits with several revolutions 
around the central BH, the most intensive GWs are radiated 
at twice the orbital frequency. It is also useful to get the in- 
verse relationship from eq. (29): 



MM,f ): 



1 



4ttGM/o 



2/3 



(30) 



According to Eq. (29), the non-relativistic result for u>o is 
adequate even in this regime. Using (13) we get 

27T(2GM) 4 / 3 _ 2/3 l-(^ /u>M,max) 2/3 ^0 
3 V °c [l-^W^max) 2 / 3 ] 

analogous to (15) for parabolic orbits, where 

2ttc 3 



dcr: 



M,max — 



32GM 



(32) 



is the maximum angular frequency, corresponding to A = 2. 

It is desirable to calculate the partial event rates of PEs with 
minimum separations bo exceeding XRsh, we substitute in 
Eq. (13), and impose the resulting constraint on the charac- 
teristic frequency: 



fo < /m,a = 27T 



2GM 



-A" 3 / 2 . 



(33) 



For marginally plunging orbits, A = 2, we get fa,x = /M,max- 

When adding up the total event rates for a particular fo 
(Eqs. 16 and 19) only the masses satisfying the constraint (33) 
have to be included in the mass integrals. Repeating § 5.1.2 
with these modification, we get for Model I 



/ if f0 V 2/3 l-2(/ //M !m ax) 2/3 dfo 



(34) 



JmJ [l-(/ //^max) 2 / 3 ] f0 

for fo < /m co ,A, and dv 1 = otherwise. For Model II, we get 
df n = / / dmidm 2 gco(m)gco(m 2 ) x 

U, , 1 ~ 2(/ //M,max) 2,/3 ( fo \ ^ d/<) 



[1-(V/M.max) 2 / 3 ] 2 ^^00^ fo ' 

(35) 

In Eqs. (34) and (35) v\ and vf{m\ , m 2 ) are the non-relativistic 
terms given by (17) and (20). 



Fig. 3 shows the resulting total event rates for dln/b in- 
tervals. The solid and dashed lines represent the total event 
rates of PEs for Models I and II, respectively, including the 
relativistic correction for encounters that avoid collisions. 
As a comparison, dotted lines display non-relativistic re- 
sults. Compared to the non-relativistic results, event rates 
decrease for two reasons: first the gravitational focusing de- 
creases the cross sections of relativistic orbits, for 2100 Hz = 
/o.max[2m NS ,(A = 2)] > fo > /o,max[2m max ,(A = 6)] = 10Hz, 
and second, the plunging orbits with A < 2 are excluded from 
our estimate. The latter effect kicks in at /o.max(2m max ,2) > 
47 Hz where the highest mass BHs suffer head-on collisions. 
At /o,max(2ffi m i n ,2) = 570 Hz even the smallest BHs are cap- 
tured, and only the NS-NS PE event rate contributions re- 
main. The NS-NS partial event rates can be visualized for 
lower frequencies by extrapolating the total event rates shown 
in Fig. 3 between 570 Hz < fo < 2100Hz. The NS-NS event 
rates are clearly negligible compared to the total rates includ- 
ing BHs. Note, that our calculations use point masses valid 
for BHs only. For fo > 1500 Hz the minimum separation of 
1 .35M Q NSs decreases under <~ 20km, for which our approx- 
imation breaks down. 

Figure 3 is useful to visualize the total PE event rate per GC. 
However, only a fraction of these events can be detected, and 
this fraction depends on both the differential encounter event 
rates d 3 v / dm\dm 2 d\n fo and also the observable distance of 
the encounter. For the detection rates we shall make use of 
the infinitesimal encounter event rate for infinitesimal mass 
and fo bins. From eq. (35) we get 11 



3 u 



dm\dm 2 dln fo 



II, \ 1 - 2(/o//M,max) 2 ^ 3 

X (mi,m 2 )- ■ —J 

[l-(/o//M,max) 2 / 3 ] 



A 

/lOO 



-2/3 



(36) 

The total event rate for one GC depends only on 
d : s v i 'dm l dm 2 d\n f , thus in eq. (35) eq. (36) was directly inte- 
grated over fo, and the mi, and m 2 distributions. However for 
the detection rate, the observation distance of the encounter 
depends on fo, m\, and m 2 differently. Therefore the differ- 
ential detection rate has a modified parameter dependence, 
implying that the integration can be carried out only after the 
observation distance had been included in the differential rate. 



5.3. Maximum Distance of Detection 

We now derive the maximum detectable distance of an 
encounter for fixed masses m\, m 2 , and characteristic fre- 
quency fo, for a given signal-to-noise ratio S/N. The luminos- 
ity distance can be expressed with the redshifted parameters, 
m iz = (1 +z)m ; for i = {1,2} and fo z = /o/(l +z), using the an- 
gular averaged signal-to-noise ratio (2) and signal waveform 
(5): 



85 tt 2 / 3 G 5//3 M, 5/3 £ rel (A) 



d L (m lz ,m 2 M= ^-^^y^gW,/*), 

(37) 

where E le i(X)/E m (X) is the enhancement of the GW energy 
for general relativistic orbits (Gair, Kennefick, & Larson 
2005, and see § 3 above), where A = A(M,/ ) = \(M z ,f 0z ) is 
given by eq. (30), and W(/,/o z ) is a factor depending on only 

1 1 We follow to notation of Miller (2002) and Will (2004) for the definition 
of partial event rates by not including the mass distribution gco(m). The mass 
distributions enter only when integrating for the total event rates eq. (42). 
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the frequencies 



W(fJ 0z ) = 



\ 



j mi 



fol 3 F(f/foz) 
f 2 S n (f) 



df 



(38) 



d h (m u m 2 ,f ) : 



F(x) is the dimensionless, normed GW energy spectrum de- 
fined in § 3, /min and / max are the minimum and maximum 
frequencies specific for GW detectors (see § 2). Henceforth 
we shall fix S /N = 5, but other values can be roughly obtained 
by scaling the final result on detection rates with (S/N)~ 3 (as- 
suming that the number of sources increases with di 3 ^). 

Note that eq. (37) formally depends on redshifted param- 
eters. However, since both the differential encounter event 
rate (36) and the GC model mass distribution depend on 
the comoving parameters, it is useful to revert to the co- 
moving parameters in eq. (37) and get di^(mi,m2,fo). This 
can be achieved by writing m !z = (1 + z)m; for i = {1,2} and 
foz = /o/(l+z) in (37) and making this equal to the standard 
formula */l cos (z) connecting the luminosity distance and red- 
shift in a specific cosmology (e.g Eisenstein 1997, the index 
"cos" refers to the cosmological luminosity distance-redshift 
formula in order to distinguish this from the maximum dis- 
tance (37) specific for PE encounters). Now both sides depend 
on z. Numerically solving for z gives z(»Ji,OT2,/o)- Finally 
substituting the result back in c/l C os(z) gives dh(mi,m 2 ,fo)- 

This procedure is however cumbersome in practice. It be- 
comes numerically very time-consuming when computing the 
total detection rates, which includes the evaluation of integrals 
over the parameters. Therefore we make the following essen- 
tial approximations when solving for the luminosity distance 
in (37): 

using z = if zi < 0.01 

using z = H (/l/c if Zi < 0. 1 
no approximations if 0.1 <z\ < 6 
using z = 6 if S/N > 5 forz = 6 

(39) 

On the RHS of eq. (39), z\ is the first approximation of the 
redshift, which is obtained by calculating d^ from the RHS of 
eq. (37) with no redshift, and making this equal to d^ cos (z\), 
and solving for zi- We neglect cosmology for z\ < 0.01 and 
take a Hubble constant H = 70km/s/Mpc for 0.01 < zi < 
0.1. Next, whenever zi > 0. 1 , we substitute z = 6 for the RHS 
of eq. (37) and in case this is already larger than c/l C os(6), then 
we conclude that the source is observable at z = 6 for S/N > 5 
and take <f L = c/l C os(6) as the maximum distance of observa- 
tion. We do not explore detection rates at larger redshifts, 
since then the BH mass and radial distribution might not have 
relaxed to the final state. If S/N < 5 for z = 6 then we execute 
the exact procedure without approximations for zi > 0.1 

Changing to the non-redshifted variables in Eq. (37) is even 
ambiguous in some cases. For NGLISA the signal-to-noise ra- 
tio is occasionally not a decreasing but an increasing function 
of the redshift. This happens when the signal is redshifted 
to the more sensitive range of frequencies of the detector, 
and the enhancement in sensitivity is more substantial than 
the attenuation from increasing the distance. In these cases a 
certain encounter can be observed within a certain distance, 
then increasing the distance, the encounter first becomes in- 
visible (i.e. S/N < 5), then again visible (i.e. S/N > 5) within 
a second maximum distance. This phenomenon occurs for 
NGLISA for large BH masses and near-maximum characteris- 
tic frequencies. 

Fig. 4 shows the maximum distance of sources with S/N = 5 
for m\=m 2 = 40M Q BH masses. The thick curves account for 
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FIG. 4. — Maximum luminosity distance, ^Lmax. of two BHs with mi = 
m 2 = 40Mq masses undergoing a PE, and emitting GWs that are detected 
on average with S/N = 5. The x-axis shows the emitted fo characteristic 
frequency of the flyby comoving with the host GC. The curves correspond 
to InLIGO, VIRGO, AdLIGO, LISA, and NGLISA, respectively. Thin dotted 
lines show the result for Newtonian waveforms (Turner 1977), thick lines 
account for general relativistic corrections to the GW amplitudes for close 
encounters. For frequencies larger than fM,ma\ = 7 1 Hz the minimum distance 
is under A = 2 Schwarzschild radii, for which a head-on collision takes place. 
We did not impose any restrictions on A for the non-relativistic curves. All 
curves account for the redshifting z, which is shown on the right border. We 
restrict to z < 6. For different masses, scales with roughly M 5 / 3 and the 
cutoff frequency scales with M~' . Since signals are broadband, the detectors 
have a chance to observe a broad range of fo. 



the relativistic corrections with minimum separations larger 
than A = 2, the thin dotted lines represent the non-relativistic 
results with no bound on A. The maximum frequency (32) 
corresponding to marginally colliding orbits (A = 2) is 71 Hz 
and scales with M~ x for other masses. For other masses ^Lmax 
scales with M 3 ! 5 . All curves account for the cosmological 
redshifting. The plot shows, that the non-relativistic approx- 
imation is adequate for small frequencies, but it implies a lu- 
minosity distance a factor of 2 - 3 lower than the relativistic 
calculation near the maximum fo frequency. Therefore, the 
dotted lines are useful to approximately visualize c/Lmax 

for 

lower M, when the cutoff frequency shifts to higher values. 

The enclosed volume and the observable sources are given 
by eq. (8). If neglecting relativistic and cosmological effects, 
we get V oc D 3 oc M 5 . 

5.4. Detection Rates 

In the previous sections we calculated the differential event 
rates of PEs for single GCs per infinitesimal mass and fre- 
quency bins, and computed the maximum distance of their 
detection. Here we combine these results to calculate the total 
detection rate of PEs. 

For fixed mi, m 2 , and f , the rate of GW detections of the 
corresponding encounters is the observed rate for a single GC 
times the number of observable GCs. Since there is a cosmo- 
logical redshift between the source GC and the observation, 
the single-GC rate is reduced by 1 +z: 



2,, total 



d 2 V 



1 



Q2 jySingle 



dm\dm2d\nfo \+zdm\dmid\nf{) 



N sc (m u m 2 J ). (40) 



The first term is the redshifted event rate expressed with the 
comoving event rate (36) and the second is simply 



N gc (m u m 2 ,f ) = N sc [t/ Lm ax(/o ,m u m 2 )] , 



(41) 
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FIG. 5. — The expected PE detection rate per logarithmic emitted char- 
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which the units of f drops out. For large f , there is an abrupt cutoff in the 
detection rate as the larger mass BHs suffer head-on collisions, leaving only 
the contribution of low-mass BHs in the PE rate. 
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FIG. 6. — The expected PE detection rate per logarithmic total mass bins, 
dlnM, for the various detectors. Note that the mass function of our GC 
model is constant for dlnm intervals for BHs 120Mq > 2m > IOMq, and 
is Gaussian type for NSs 2m ~ 2.7Mq . BH-NS encounters dominate for 
6.35Mq < M < IOMq. The dominant PE contribution is expected from 
m\2 =40-60Mq component masses. 



where N sc (d\J is the number of GCs within a given maximum 
luminosity distance. In practice t/Lmax is given by eq. (39) 
which we substitute in eq. (8). 

The total detection rate is then simply the integral of the 
differential detection rate (40) using the CO mass distribution. 
After substituting, we get 



rfo« 
Jo 



fo 



fo<fn 



gco(mi)gco(m2)- 



Am\Ani2 



<9V< 



(42) 



dm\dni2d\n fo 

where gco(m) is the CO mass distribution for Model II, de- 
fined in § 4.3.2. The mass integrals are evaluated over the 
(mi,m2) domain for which the encounter avoids a collision 
(i- e - fo < fM,mnx, see eq. [32]) and the fo integral extends 
to a maximum possible frequency independent of masses 
(~ /2m NSl max)- The result of eq. (42) is one number, the ex- 
pected rate of detection for the specific detector. 

5.5. Results 

The estimated total number of successful detections from 
eq. (42) is v m = 5.5 x lO^yr" 1 for InLIGO, 7.2 x lO^yr" 1 
for VIRGO, 6.3 x lO^yr" 1 for AdLIGO, 2.9 x lO^yr" 1 for 
LISA and 1 .Oyr" 1 for NGLISA. 

It is interesting to see the differential event rate per logarith- 
mic fo bin independent of masses, which is obtained by carry- 
ing out only the mass integrals in eq. (42). The result is shown 
in Fig. 5. The figure shows that both AdLIGO and NGLISA 
could have some chance to detect PE events, if observing 
for one year; AdLIGO mainly sensitive to /o frequencies be- 
tween 30 and 80 Hz, and NGLISA sensitive between 0.2 Hz 
and 10 Hz. There is a sharp cutoff in the PE detection rate 
for high frequencies. In this regime, the encounters among 
the relatively higher mass BHs are not parabolic, but result 
in direct captures, and only the lower mass BHs contribute 
to the PE detection rate, show the differential event rate for 
logarithmic total mass bins, dlnM, and for logarithmic mass 
ratio bins, din q. Here we define the mass ratio as q = m < /m > 
for which q < 1. (Recall the definitions m< = min(mi,m2) 



and m> = max(mi,m2).) The M dependent partial PE detec- 
tion rate can be obtained from eq. (42) by changing the mi,ni2 
integrals to M and ni2 variables, rearranging the order of inte- 
grals, and evaluating the fo and ni2 integrals only. The partial 
PE rates for fixed q can be obtained similarly, by changing 
to m> and q variables, and evaluating the /o and m> inte- 
grals only. In § 5.1.2 we demonstrated that the event rates 
of GCs are sensitive to (m 10 / 3 ), and are inclined towards the 
high-mass end of the CO distribution, in particular PEs of NSs 
have a relatively negligible event rate. The detectable volume 
entails an even stronger mass dependence m 5 . Therefore, for a 
mass distribution of m" 1 , we expect a scaling with <~ m 22 / 3 for 
logarithmic total mass bins, implying that the highest mass 
BHs will dominate the PE event rates. However, increasing 
the BH masses decreases the maximum fo frequency of an 
encounter avoiding a collision. Figure 6 verifies that all of 
the detectors are indeed much more sensitive to large total 
masses, even though our model GC (i.e. Model II) contains a 
small relative number of these objects. Note, that the BH mass 
distribution ^bh(^) is constant for dlnm bins. With Figure 6 
the partial detection rates of BH-BH and NS-NS encounters 
can be visualized. For M > 2m min = 10M Q the BH-BH en- 
counters dominate, while M « 2mNs correspond exclusively 
to NS-NS encounters. PE detections of NS-NS encounters 
are practically impossible, they are suppressed by at least 9 
orders of magnitudes. Similarly, Figure 7 shows that BH-NS 
encounters are also suppressed by 6 orders of magnitudes! 

Figure 8 shows the detection rate as a function of mini- 
mum distance, Ami n of the encounters. Recall that for a given 
total mass M, A determines the characteristic frequency fo 
by eq. (33), and marginally plunging orbits correspond to 
A m ; n = 2. Figure 8 was obtained by changing the domain of 
integration of fo to fo < /m,a m ec l- (42). The curves show that 
terrestrial detectors are more sensitive to close approaches 
than space detectors. The A m ; n = 2 case corresponds to all of 
the PE detections. It is interesting to note, that terrestrial de- 
tectors display a different A dependence: AdLIGO rates show 
a weaker increase for marginally colliding orbits A ~ 2. This 
is a consequence of cosmology: the observation distance is so 
large (Fig. 4) that the cosmological comoving volume element 
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Fig. 7. — The expected PE detection rate as a function of the mass ratio 
q — m^lmx of the interacting masses, where m\ > ni2 is assumed. The par- 
tial PE detections are plotted per logarithmic mass ratio bins, din*?, for the 
various detectors. Note that the assumed smallest and largest CO masses of 
our GC model implies a cutoff below 1.35Mq /6OM0. The detection rate is 
dominated by equal mass encounters. PEs with q > 5M0/6OM0 are dom- 
inated by BH-BH encounters, while 1.35Mq/60Mq < q < 5M /6OM 
correspond to BH-NS events. 
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FIG. 8. — The total cumulative expected detection rate of PEs subject 
to the constraint that the minimum separation exceeds A times the total 
Schwarzschild radius. For A < 2 a head-on collision occurs, which we are 
not considering in the present paper (see § 7.1 for a discussion). 



is significantly smaller, and the GW frequency is redshifted 
outside the sensitive domain of the detector for these events. 

6. CONCLUSIONS 

PEs of solar BHs are possible sources of gravitational ra- 
diation. Our results suggest that current and near future GW 
detectors are potentially capable of detecting these signals in 
the local universe and up to cosmological distances for the 
higher masses. We anticipate S/N > 5 matched-filtering de- 
tection rates for quasi-parabolic trajectories avoiding colli- 
sions, v lot = 5.5 x lO^yr" 1 for InLIGO, 7.2 x lO^yr" 1 for 
VIRGO, 0.063 yr" 1 for AdLIGO, 2.9 x lO^yr" 1 for LISA, and 
l.Oyr" 1 for NGLISA. For different signal-to-noise ratios, de- 
tection rates scale by approximately (S/N)~ 3 . These results 
correspond to a BH mass function gBH(m) oc m" 1 with mini- 
mum and maximum masses of 5 and 60M Q . For comparison 
we ran calculations for more general distributions ^bh(wi) oc 
m~ p , and different m m \ n and m max values. Table 1 lists the re- 
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[M ] 
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[yr '] 
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[yr '] 


r -In 

[yr '] 
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[yr '] 
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i o r— 6i 


5 2 f— 71 
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60 


5 


2 2 (-41 


3 1 C— 41 


2 4 (-11 


1 2 (-51 


5.2 


n 


60 


40 


i.jy—**) 


2 1 C 41 


1 3 (-11 


7 9 (-61 


4.3 
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5 


i.vy j) 


3 .6(-3) 


1.3 


1 3 (-41 
i .j y^+f 




Q 


100 


40 


9 5 T — 41 


3 4 f— 31 


1.2 


1 2 (-41 
± .z. y^tf 




Q 


100 


80 


2 4 (-41 


1 2 C— 31 


3 4 (-11 


4 3 (-51 




[ 


20 


5 


4.8 (-7) 


2.4 (-7) 


1.0 (-3) 


4.5 (-8) 


2.7 (-4) 




60 


5 


5.5 (-5) 


7.2 (-5) 


6.3 (-2) 


2.9 (-6) 


1.0 




60 


40 


2.5 (-5) 


3.9 (-5) 


2.6 (-2) 


1.3 (-6) 


7.6(-l) 




100 


5 


2.2 (-4) 


6.4 (-4) 


2.8 (-1) 


2.4 (-5) 






100 


40 


1.8 (-4) 


5.7 (-4) 


2.2(-l) 


2.0 (-5) 






100 


80 


2.9 (-5) 


1.4 (-4) 


4.2 (-2) 


5.2(-6) 




2 


20 


5 


1.8 (-7) 


9.1 (-8) 


4.1 (-4) 


2.1 (-8) 


9.5 (-5) 


2 


60 


5 


6.6 (-6) 


7.2 (-6) 


8.3 (-3) 


3.7 (-7) 


8.4 (-2) 


2 


60 


40 


1.8 (-6) 


2.8 (-6) 


1.9 (-3) 


1.0 (-7) 


5.0(-2) 


2 


100 


5 


1.7 (-5) 


3.6 (-5) 


2.1 (-2) 


1.5 (-6) 




2 


100 


40 


9.7 (-6) 


2.7 (-5) 


1.1 (-2) 


9.6 (-7) 




2 


100 


80 


8.9 (-7) 


4.3 (-6) 


1.3 (-3) 


1.6 (-7) 





NOTE. — The number of BHs per GC is normalized to Wbh = 500 for 
m min = 5Mq for all choices of m max . For larger m min , all BHs with 5M < 
m < m mm are assumed to have escaped from the cluster. Detection rates are 
given in normal form, where the exponent is shown in parenthesis. Some 
fields left blank correspond to cases where the S/N is not a monotonically 
decreasing function of distance (see § 5.3). 



suits for these models. Here, we have fixed A^bh = 500 for 
models with m m ; n = 5M Q . For larger m m ; n , we reduce Nbh 
by assuming that BHs with masses 5M0 < m < m m ; n have 
escaped the cluster. Results are very different for various 
choices of parameters (see § 7.3.2 for a detailed discussion 
below). 

We constructed two different GC models. We conclude that 
a uniform mass and density distribution (Model I) is inade- 
quate since the contribution of the GC core consisting of the 
more massive BHs are significantly underestimated. After ac- 
counting for mass distribution and mass segregation, as well 
as the relative velocity distribution of the sources (Model II) 
we obtained event rates two orders of magnitudes higher than 
Model I per GC. Moreover, more massive BHs in GCs are 
visible to significantly larger distances, and supply the most 
prominent sources of PEs for detection (Fig. 6). 

In § 7, we include a critical review of our assumptions and 
their influence on the results. To point out just one thing, note 
that compared to our previous estimates above, the PE detec- 
tion rates might have been underestimated by four orders of 
magnitudes for core-collapsed GCs (depending on the final 
core radius and population, see § 7.3.3 below)! 

7. DISCUSSION 

7.1. Comparison with Other Orbits 

How do PE event rates compare to the event rates of binary 
inspirals, mergers, and head-on collisions? What are the main 
factors for the difference? We briefly discuss these questions 
in this section. 

7. 1. 1. Basic Features of Parabolic Encounters 

Let us quickly summarize the main properties of PE 
sources. 

• The event rates for din fo intervals scale with f^ 2 ^ for 
trajectories avoiding collisions. Collisions decrease PE 
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event rates quickly for large frequencies, / > 50 Hz, 
rates drop four orders of magnitudes between 50 and 
500 Hz. 

—1 /3 

• The /o-scaling of the amplitude of the signal is / ' 
for GW frequencies / <~ fo. The integrated RSS signal 

amplitude scales with f^ b . 

• The signal energy-spectrum is broadband, has a max- 
imum at / <~ fo, and a relatively shallow cutoff for 
larger frequencies. The half-maximum of the signal is 
Af <~ 1.5/o, and the spectral energy density drops to 
1% at /~ 5/ (Turner 1977). 

• In terms of detections for dln/b intervals, the maxi- 
mum distance of PE detections for a band-pass detector 
is roughly independent of frequencies between / m ; n < 
fo < /omai, where / m i n is the minimum detectable fre- 
quency of the detector and /o max is the maximum char- 
acteristic frequency of a PE avoiding collisions. 

• The detection rates of equal-mass PEs scale with 
(m 22 / 3 ) for dlnm intervals, the higher mass BHs domi- 
nate PE detections. 

• Space detectors will possibly detect more PE events in 
the local universe, but terrestrial detectors see further. 

• Typical event rates are 1.6 x 10" 12 yr _1 GC _1 or equiva- 
lent^ TZ = 1.4 x lO-'^yr" 1 Mpc -3 (A In /(i)" 1 for f = 
50 Hz. 

• Typical maximum distance of detection for PEs with 
appropriate f is <~ 300 Mpc for InLIGO and VIRGO, 
z - 1 for AdLIGO, - 0.4 Mpc for LISA, and z ~ 0.2 for 
NGLISA. 

• Typical overall rate of PE detections per year is ~ 10" 4 
for InLIGO and VIRGO, - 0.1 for AdLIGO, - 10" 6 
for LISA, and - 1 for NGLISA. 

Table 2 contrasts the event rates of PEs and other possibly 
detectable sources in the literature. The event rates, TZ, listed 
in the table are normalized to an average space density of 
galaxies n ga i = 0.029 Mpc -3 and n gc = 2.9 Mpc -3 for GCs (see 
§ 4.2). For PEs the TZ values shown correspond to the sensi- 
tive range of characteristic frequencies: 10Hz < fo < 200 Hz 
for InLIGO, VIRGO, AdLIGO, 10" 4 Hz < /„ < 1 Hz for LISA, 
and 10" 1 Hz < f < 60 Hz for NGLISA (see Fig. 5). The max- 
imum distance of detection, d^, shown are typical values for 
detectable encounters. The elements marked with "-" are not 
given by the corresponding references. Values are left blank 
that are the same as in the previous row. Note that the num- 
bers are very uncertain depending on model assumptions. PE 
event rates correspond to our standard GC model, most signif- 
icant uncertainties are the numbers of higher mass BHs in the 
cluster. Note, that stellar BH-BH inspiral event rate estimates 
in GCs vary 3 orders of magnitudes! 

7.1.2. Parabolic Encounters vs. Inspirals 

The event rates of PEs depend on the characteristic fre- 
quency. Without any specifications PEs are more regular than 
e.g. BH-BH inspirals. However, within the sensitive range 
of GW terrestrial detectors PEs are rather rare by a factor of 
~ 10" 4 . Event rates are higher for space detector frequencies, 



TABLE 2 

Comparison of event and detection rates with other sources 



Event 


Loc 


'Tl 

f-C 


Detector 


j 

"L 

[Gpc] 


V 

[yr 1 ] 


BH/BH PE 1 


GC 


7 H-ll) 


InLIGO 


3 0(-ll 


5.5 (-5) 




GC 


7.1 (-1 1) 


VIRGO 


3.0(-l) 


7.2 (-5) 




GC 


7.1 (-1 1) 


AdLIGO 


5.6 


6.3 (-2) 




GC 


2.4 (-7) 


LISA 


5.1 (-4) 


2.9 (-6) 




GC 


2.4 (-9) 


NGLISA 


1.9 


1.0 


NS/NS insp 2 


Field 


1.9 (-7) 


AdLIGO 


2.0(-l) 


3.0 


( 3 LPP) 


Field 


1.1 (-6) 


AdLIGO 


2.0(-l) 


3.4(+l) 


( 4 Kalogera) 


Field 


2.0 (-6) 


InLIGO 


2.0 (-2) 


3.5 (-2) 






2.0 (-6) 


AdLIGO 


3.5 (-1) 


1.9 (+2) 


( 5 NGF) 


Field 


2.9 (-4) 


InLIGO 


2.0 (-2) 


3.4(-l) 






2.9 (-4) 


AdLIGO 


3.0(-l) 


1.1 (+4) 


BH/NS insp 5 


Field 


2.9 (-4) 


InLIGO 


4.3 (-2) 


3.3 (+1) 






2.9 (-4) 


AdLIGO 


6.5 (-1) 


1.2 (+5) 


BH/BH insp 6 


Field 


2.0 (-9) 


AdLIGO 


1.1 


1.0 


( 7 PZM) 


Nucleus 


1.8 (-8) 


AdLIGO 


1.1 


1.0 (+2) 




Zero-age GC 


3.2 (-7) 


AdLIGO 


1.1 


1.8 (+3) 




Evolved GC 


5.4(-8) 


AdLIGO 


1.1 


3.0 (+2) 


( 8 Miller) 


GC 


3.5 (-9) 


AdLIGO 


1.6 


1.0(+1) 


('O'Leary) 


Zero-age GC 


8.7 (-7) 


AdLIGO 


- 


2.7 (+3) 




Evolved GC 


8.7 (-10) 


AdLIGO 




2.7 


IMBH/BH 8 


GC 


3.5 (-9) 


AdLIGO 


1.6 


4 0(+l) 






3.5 (-9) 


LISA 


2.0 (-1) 


7.0 (-3) 


( 10 Will) 


GC 


8.5 (-9) 


LISA 


4.0 (-2) 


1.0 (-6) 


( n GMH) 


GC 


5.1 (-5) 


LISA 


4.0(-2) 


6.0 (-3) 


( 12 HPZ) 


Field 


1.4 (-7) 


LISA 


3.5 (-1) 


8.6 


SMBHAVD 13 


Nucleus 


1.7 (-8) 


LISA 


1.6 


9.4 (+1) 


( 14 HA) 


Nucleus 


9.0 (-8) 


LISA 


1.0 


1.3 (+2) 


SMBH/BH 13 


Nucleus 


3.2 (-9) 


LISA 


6.6 


1.1 (+3) 


SMBH/IMBH 13 


Nucleus 


2.9(-12) 


LISA 


6.6 


1.0 


( 15 PZ) 


Nucleus 


8.3 (-10) 


LISA 


4.4 


1.0 (+2) 


NOTE. — Event rates TZ are is in units ft 3 yr 


4 Mpc~ 3 . 







1 Parabolic encounters for our standard GC model. The maximum distance of 
detection corresponds to mi = ni2 = 5OM0, fo = 50 Hz for InLIGO, VIRGO, 
AdLIGO and NGLISA, and f = lOmHz for LISA. 

2 NS/NS inspirals, Portegies Zwart & Spreeuw (1996) theoretical models of bi- 
nary evolution calibrated to the observed supernova rate. 

3 Lipunov, Postnov, & Prokhorov (1997) 

4 Kalogera et al. (2004), based on 3 highly relativistic radio pulsars. 

5 Nakar, Gal- Yam, & Fox (2005), based on recent observations of 4 short/hard 
gamma-ray bursts, identified with NS/NS or BH/NS inspirals. 

6 BH/BH inspirals, Portegies Zwart & Yungelson (1998). 

7 Portegies Zwart & McMillan (2000). Distance corresponds to the inspiral of 
masses m\=ni2 = 10M Q and detection rates correspond to a 100% detection effi- 
ciency within this distance. 

8 Miller (2002) 

9 O'Leary et al. (2006). They compute the number of mergers as a function of 
time, as BHs are being slowly depleted from the cluster. The detection rates cor- 
respond to cases, where all GCs are zero-age (lMyr) or evolved (lOGyr), with a 
number density 1 Mpc~ 3 . 

10 Stellar BH inspirals into intermediate mass BHs, adopted from Will (2004). 
Distance corresponds to the inspiral of masses m\ = IOOMq and mi = IOMq, TZ 
corresponds to the inspirals with a maximum time 40^100 yr before merger which 
can be detected with S/N =10 after a 1 yr integration with LISA. 

1 1 Giiltekin, Miller, & Hamilton (2006) assuming more optimistic probability of 
IMBHs in GCs and IMBH maximum mass than Will (2004) 

12 Hopman & Portegies Zwart (2005), assuming observed ultraluminous X-ray 
sources are associated to tidal captures of IOMq BHs by 10 3 Mq IMBHs. 

13 Inspirals into supermassive BHs, M = 10 6 Mq . Gair et al. (2004) assuming their 
optimistic set of results. 

14 Hopman & Alexander (2005) 

15 Portegies Zwart et al. (2006) N-body simulations of the formation and inward 
migration of IMBHs in the Galactic center. 
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however space detectors have a smaller distance of maximum 
observation. 

Table 2 suggests that the detectable distance of PEs is com- 
parable to inspirals. However, this is somewhat deceptive 
since PE results correspond to larger BH masses, 50M Q rather 
than 1OM which is regular for BH-BH inspirals in the liter- 
ature. For IOMq component masses, the maximum distance 
of observation is less for PEs than for inspirals (Fig. 4), the 
PE event rates are suppressed by a factor of ~ 10 3 (Fig. 6). 
However for larger masses the comparison changes with the 
following factors. First, although the GW signal amplitude is 
proportional to A4 5 ^ 3 for inspirals, increasing the masses re- 
duces the signal frequency (which in turn reduces the detector 
sensitivity) and also reduces the observation time (which also 
decreases the effective signal amplitude). Another important 
difference is in the binary separation A which determines the 
signal frequency. It is restricted to A > 3 for inspirals, the in- 
nermost stable circular orbit, a more stringent constraint than 
the condition A > 2 for PEs (see § 5.2). Finally opposed to 
the PE signal waveforms, the inspiral signals are narrow band, 
implying that the high mass, i.e. low frequency, inspiral wave- 
forms are much harder to detect as much smaller signal power 
accumulates at the more sensitive range of frequencies. When 
combining all of these effects we expect that low mass BH in- 
spirals are detectable further with terrestrial detectors, while 
for large masses where the observation is limited to at most 
a few orbits, marginally collisional PEs are better detected. 
Therefore, detectors can observe the higher mass encounters 
for PEs. This is exactly analogous to the comparison of the 
inspiral and plunge phases of binary coalescence, for which 
the detection of plunge dominates for large masses (Flanagan 
& Hughes 1998). 

Among the GW detection candidate sources in GCs, PEs 
are very infrequent compared to stellar BH-BH inspiral rate 
estimates of Portegies Zwart & McMillan (2000) or Miller 
(2002) within GCs, but are comparable to the recent results 
of O'Leary et al. (2006). Observations of radio pulsars and 
gamma ray bursts suggest several orders of magnitude larger 
numbers for NS-NS or BH-NS inspiral detections (Kalogera 
et al. 2004; Nakar, Gal- Yam, & Fox 2005). 

7.1.3. Parabolic Encounters vs. Head-on Collisions 

As a second example, let us consider the event rates of 
head-on collisions for unbound encounters. Head-on colli- 
sions are related to PEs, by extending the parameter A to val- 
ues less than 2, the unstable circular orbit. Thus it would 
be relatively straightforward to extend the analysis to these 
events, by examining the event rates for small initial impact 
parameters, and computing the detectability as a function of 
this parameter. However, the exact shapes of gravity wave- 
forms are presently not available for collisions (see Baker et 
al. 2006 for current progress), therefore maximum likelihood 
detections are not possible, and the detection of these bursts 
requires much higher signal-to-noise levels. We shall argue 
that direct head-on collision detections are potentially less fre- 
quent than PEs. 

The rate of head-on collisions between BHs is well known 
(e.g Hills & Day 1976; Cutler, Kennefick, & Poisson 1994; 
Sigurdsson & Rees 1997) however the detection rates of the 
resultant GW signals is subject to the uncertainty of the GW 
signals (Flanagan & Hughes 1998). Direct collisions produce 
potentially less intensive GW signals than close PEs even if 
neglecting the relativistic amplitude enhancement for PEs. To 
see this, let us compare GW signal strengths that we adopt 



for PEs (see § 3, and Turner 1977; Martel 2004; Gair, Ken- 
nefick, & Larson 2005), with general relativistic calculations 
for head-on collisions. To our best knowledge, off-axis col- 
lisions of BHs have not been calculated as a function of im- 
pact parameter. For radial head-on collision of Schwarzschild 
BHs AE ~ 0.01 (fi 2 /M)c 2 (Davis et al. 1971; Anninos et al. 
1993; Moreschi 1999; Sperhake et al. 2005). Sasaki & Naka- 
mura (1982) derived GW energies for the radial infall of a test 
particle into a Kerr BH, and Mine, Shibata, & Tanaka (1996) 
accounted for the spin of the infalling particle in addition. Re- 
sults are in the range of AE = 0.03 - 0.01 (/x 2 /M)c 2 accord- 
ing to the magnitude and alignment of spins and the relative 
direction of the approach. For high-velocity head-on colli- 
sions, there are significantly larger results: AE = 0.328 /ic 2 
for non-rotating BHs (D'Eath & Payne 1992), and up to 
AE = 0.70 ^c 2 for extreme Kerr-BHs (Cardoso & Lemos 
2003). However, in GCs the initial velocities are typically 
non-relativistic, therefore we do not expect a significant rel- 
ative contribution of relativistic head-on collisions. Once the 
BHs are so near that a common surrounding horizon envelope 
forms the space-time relaxes to a Kerr-BH. The energy out- 
put of this process is between AE = 6 x 10~ 6 Mc 2 (Price & 
Pullin 1994) for axisymmetric encounters and < O.OlMc 2 for 
quasi-circular initial conditions (Khanna et al. 1999). In com- 
parison, the energy output in GWs for non-relativistic PEs is 
A£ PE = 0.01 (A/4.1)" 7 / 2 (^ 2 /M)c 2 (Turner 1977). Using the 
low-velocity case, the GW amplitudes of BH collisions are 
overestimated by the Newtonian results by a factor between 
(A/3)" 7 / 2 and (A/4)" 7 / 2 , depending on spins. Therefore, the 
extapolation of the Newtonian treatment to the regime where 
the minimum separation is A <C 3 leads to significant overes- 
timates of the true head-on collision GW energies. In con- 
clusion, the extrapolation of event rates as a function of A (for 
A < 2 in Fig. 8) or as a function of the logarithmic characteris- 
tic frequency (for /o > f\=2 in Fig. 4, dotted lines) is possibly 
overly optimistic and therefore inadequate for the estimation 
of the detection rates of head-on collisions with the particular 
GW detectors. 

7.2. Approximations in the Analysis 

Our event rate estimates rely on several approximations. 
The most important caveat in our analysis is possibly neglect- 
ing GW recoil capture in bound eccentric orbits. The GW 
radiation reaction is substantial for strong gravitational fields, 
for low A. For initially nearly parabolic orbits, the perias- 
tron distance and the eccentricity is decreased 12 (Cutler, Ken- 
nefick, & Poisson 1994). The first consequence is a minor 
decrease in the PE event rate, because of the increase of the 
cross-section of direct capture. On the other hand, GW re- 
coil produces bound orbits from initially unbound trajectories 
(Lee 1993). The periastron distance is then further decreased 
during each subsequent close approach inducing successively 
stronger GW radiation. Therefore PE events are potential pre- 
cursors of multiple subsequent more intense highly eccentric 
bound encounters, analogous to the captures of stellar com- 
pact objects by supermassive black holes (Hopman & Alexan- 
der 2005). The GW detection rate of the resultant orbits is 
likely to be significantly higher than PE detections. As a re- 
sult, we anticipate several successful detections for AdLIGO 
per year for a wide range of BH mass-distribution models (see 

12 unless A > 2.05, in which case the eccentricity is increased by GW 
recoil 
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Tab. 1 for PEs without GW recoil capture). We leave a de- 
tailed quantitative study for a future paper. 

There is a second independent reason suggesting that higher 
detection rates will be more likely. Throughout the paper, we 
estimated matched filtering detection signal-to-noise ampli- 
tudes with the angular averaged formula which is valid for an 
analysis using only a single GW detector. However, a coin- 
cident analysis incorporating several detectors allows much 
more optimistic detection limits (see Jaronowski et al. 1996, 
and § 3 above). If lowering the angular averaged minimum 
detection limit to S/N = 3/V5 (equivalent to an optimal- 
orientation single-detector observation at S/N = 3) yields 
i/ ot = 1.7 x 10 _3 yr _1 forlnLIGO, 2.7 x lO^yr" 1 for VIRGO, 
0.46 yr" 1 for AdLIGO, and 9.5 x lO^yr" 1 for LISA for our 
standard GC Model II with N BU = 500, p = 1, m min = 5M Q , 
and m max = 60M Q . For models with larger m max detection 
rates are even higher (see Tab. 1), implying several successful 
detections per year for advanced terrestrial detectors. 

We have also neglected the bound binary interactions in 
the scattering dynamics and restricted only to single-single 
encounters. Depending on the angle of injection this could 
increase or decrease event rates. However our analysis of 
single-single interaction shows (§ 5.1.2), that the typical nu- 
merical values for PE cross-sections are extremely small, so 
that the injection has to have an initial velocity very accu- 
rately pointed towards the target CO in order to produce a de- 
tectable signal. For typical encounters with m\=m-i = 50M Q 
and /o = 50 Hz, the minimum separation is bo <~ 10~ 6 AU, 
and the impact parameter IS boo ^ 

10 _2 AU. As a result, for 
single-binary interactions we speculate that the separation of 
scales is possible to distinguish three independent phases of 
the interaction: (i) the faraway zone r a^, (ii) the inter- 
mediate zone r ~ flbi n , and (iii) the PE zone r ~ bo <C flbin- 
In (iii) the binary companion can be discarded. Moreover, 
note that the velocity during (ii) is still negligible compared 
to (iii). Therefore in practice, the beginning of phase (iii) is 
exactly analogous to the initial conditions of a single-single 
encounter. The only difference is the distribution of velocities 
is not isotropic, but after phase (i) it is beamed toward the cen- 
ter of mass, and phase (ii) adds a random deflection due to the 
companion. Plugging in the numbers for binary separations of 
flbin 3> 10 AU we conclude, binary focusing is not likely to 
significantly modify our PE rate estimates. Numerical simu- 
lations would be needed to determine the exact modifications 
in the estimates. 

Binary interactions also alter the total number and mass dis- 
tribution of BHs in the cluster. However in our calculations 
the total number and mass function of BHs are input parame- 
ters, which can be chosen consistently with the most sophisti- 
cated simulations. 

Throughout our analysis we assumed simplified GC mod- 
els. While our most sophisticated model accounts for the 
mass distribution, mass segregation, and relative velocities 
(see § 4.3) it does not consider the nonuniform radial distri- 
bution of density of regular stars in the cluster core, nor does 
it consider variations around the characteristic GC model pa- 
rameters (e.g. virial radius, total mass, etc). However, the 
final results are simple powers of the characteristic parame- 
ters (v oc q 2 NtjR~ 2 5 ). Our treatment allows upper and lower 
bounds to be made on the exact GC model detection rates. 
These bounds are still much tighter than other sources of un- 
certainties, which justifies the simplifying model assumptions 
in this analysis. 



Another major approximation was to adopt the angular av- 
eraged signal waveforms in the Newtonian approximation 
(Turner 1977), and corrected for the relativistic enhance- 
ment of the amplitude, substantial for close-encounters. We 
adopted the relativistic correction for the quadrupole radia- 
tion of a test particle geodesies (Martel 2004; Gair, Kennefick, 
& Larson 2005) and extrapolated results for other masses. 
These estimates do not account for GW recoil. However, 
Fig. 8 shows that the contribution of extreme zoom-whirl or- 
bits A w 2 does not ruin our estimates, since the detection rate 
does not increase substantially for marginally plunging orbits. 
GW recoil reduces PE signal-power by driving the interacting 
masses to collisions, thereby terminating extreme zoom-whirl 
orbits much sooner than the no-recoil encounter time (Gair, 
Kennefick, & Larson 2006). We conclude that neglecting GW 
recoil did not lead to a large overestimate, implying that our 
results are acceptable approximations in this respect. 

An exact treatment would have to utilize the more exact 
post-newtonian waveforms of the general problem using ar- 
bitrary masses and spins, and should take into account the 
forward peaking of GWs for high velocities, Doppler shift of 
GW frequencies, spin-orbit, and spin-spin interactions, etc. 
Although it is clear that a real data analysis matched filter- 
ing would have to be carried out with exact signal templates, 
the leading-order (i.e. Newtonian) term dominates the angu- 
lar averaged signal power, which is therefore an adequate first 
estimate for the detection rates. 

7.3. Uncertainties in the Result 

7.3.1. Model Parameters 

There are several uncertainties in our estimate. Among 
the most important uncertainties are the values of the GC 
model parameters, like the number of BHs in the cluster A^h- 
Portegies Zwart & McMillan (2000) derives q Bli = A^h/AW = 
6 x 10 -4 by using Scalo (1986) initial mass function (IMF) 
and assumed that every object more massive than 20M Q up 
to 100M Q had evolved to a BH. When using a Salpeter IMF, 
the result is q Bii = 10" 3 (Miller 2002), and Kroupa & Weid- 
ner (2003) IMF gives q BU = 1.5 x 10" 3 (O'Leary et al. 2006). 
We adopt the most conservative result of Portegies Zwart & 
McMillan (2000). However there is a chance that a non- 
negligible fraction of the stars have been ejected from the 
cluster or have undergone subsequent mergers. Both pro- 
cesses increase the estimate on the final BH fraction (Miller 
2002). On the other hand dynamical binary interactions, bi- 
nary recoil kicks, or GW recoil of BH mergers can eject 
BHs, thereby reducing their overall numbers and possibly also 
modify the mass-distribution. In fact, a significant portion of 
the stellar-mass BH population might be ejected, especially in 
small clusters (Sigurdsson & Hernquist 1993; Portegies Zwart 
& McMillan 2000; O'Leary et al. 2006). Belczynski et al. 
(2005) find that for an initial binary fraction of 50%, the re- 
tained fraction of BHs varies between 0.4 and 0.7. In our fidu- 
cial calculations we adopted q BH = 5 x 10" 4 and N tot = 10 6 . To 
see the effects of BH ejection, Tab. 1 shows results for other 
models. For the general case, we provide analytical scalings 
which can be readily used in case these parameters are bet- 
ter determined in the future. For example since v lot oc A^ H , 
for A^bh = 5000 (50) detection rates increase (decrease) by a 
factor of 100. 

7.3.2. Black Hole Mass Distribution 
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An even more significant source of uncertainty is the mass 
distribution of BHs in the clusters. We have calculated detec- 
tion rates for several distributions (Tab. 1). Increasing the p 
exponent of the distribution, g BH oc m~ p , decreases the detec- 
tion rate by a factor of ~ 5 for a unit change in p. Chang- 
ing the maximum mass of the distribution varies the results 
even more significantly. Compared to the detection rate corre- 
sponding to m max = 6OM , for m max = 100M Q (20M Q ) we get 
a <~ 1 - 3 order of magnitude increase (decrease) depending 
on the detectors. We have also tried changing the minimum 
mass m m ; n , by assuming that the BHs with masses m < m m ; n 
have escaped from the cluster. Compared to m m ; n = 5M Q , a 
value of m m ; n = 40M Q reduces detection rates by a factor of 
~ 2. In the appendix, we provide ready-to-use formulas for 
calculating detection rates for other parameter values. 

From the theoretical point of view, simulations of the initial 
stellar BH mass function (Fryer & Kalogera 2001) result in a 
maximum mass limit of ~ 20M Q , but the particular form of 
the mass function is very different for various assumptions 
(fraction of explosion energy used to unbind the star, stel- 
lar winds, mass transfer after helium ignition, etc). Recent 
simulations of rapid star evolution assuming a lower metal- 
licity for the progenitor stars (weaker stellar winds) appropri- 
ate for GCs and including a large fraction of binaries, colli- 
sions, and accretion leading to the mass buildup of BHs im- 
ply a stellar-mass BH distribution with maximum BH masses 
around M max = 60 - 100M Q (Belczynski et al. 2005). Simu- 
lations of the subsequent long-term dynamical evolution has 
been shown to be sensitive to BH binary and triple interac- 
tions (O'Leary et al. 2006). Binary-single body interactions, 
BH-star collisions, and GW recoil kick can possibly signif- 
icantly reduce the low mass BH population but enhance the 
mass of the most massive BHs in the cluster. From the obser- 
vational point of view, there is yet lacking evidence for stellar 
mass BHs with m > 20M Q , but this might be accounted for 
the low number statistics (a total of 20 X-ray stellar-mass BH 
candidates have been identified to date, Casares 2005). 

7.3.3. Core Collapse 

Finally, a considerable uncertainty in the PE detection rates 
results from the actual scaling of the mass segregation rela- 
tionships. Even in our complicated model we have assumed 
a simple mass segregation, based on thermal equipartition 
among the specific CO components. This assumption is in 
fact valid only among the decoupled high mass components 
within the core. Spitzer (1969) has shown that in simple 
two-component systems consisting of masses m\ and m.2, with 
m\ -C tn.2, global equipartition cannot be attained if the low- 
mass component determines the potential everywhere in the 
cluster. In this case, the high mass components become dy- 
namically decoupled from the rest of the cluster, and the clus- 
ter core collapses to a much smaller radius, R cole . This pic- 
ture has been confirmed by numerical simulations for more 
general mass functions (Watters, Joshi, & Rasio 2000, and 
references therein). Giirkan, Freitag, & Rasio (2004) showed 
that approximate local thermal equipartition is attained within 
the core, and velocities follow v m = {Km/ m B \i)~ l l 2 v COK , where 
v core is the velocity dispersion, bibh is the mass of compo- 
nents in the core, K describes the departure from equiparti- 
tion, it is a number of order 1. The total time of the col- 
lapse and the final magnitude of core velocities or core radius, 
depends sensitively on the initial fraction of binaries. For a 
single mass cluster Heggie, Trenti, & Hut (2006) found that 
0.01 <R C0K /R^ < 0.1, larger values valid for a large frac- 



tion of binaries (here R gc is the half-mass radius). In contrast 
our simple mass segregation led to /?5om q = 0.l4R gc , which is 
a factor 1.4-14 higher. Note, that the virial theorem implies 
v m oc R m for a homogeneous mass distribution. Detection rates 
scale with R^v~^, and the contribution of m <~ 50M Q domi- 
nated the final results (see Figs 6 and 7). Therefore post-core 
collapse mass segregation implies detection rates increased 
by (1.4) 4 -(14) 4 . Thus, in the most optimistic case, we get 
a substantial increase in the detection rates, i.e. 2.1 yr" 1 for 
InLIGO, 2.8 yr" 1 for VIRGO, 6.6day _1 for AdLIGO, 0. 1 yr" 1 
forL/&4, and 4.4 hr" 1 for NGLISA ! 

On the other hand, if core collapse leads to runaway col- 
lisions and the buildup of a single intermediate mass black 
hole, while stellar mass BHs are ejected from the cluster (Fre- 
itag, Giirkan, & Rasio 2006), PE detection rates might be 
considerably suppressed after collapse (i.e. v oc A^ H ). More 
information on the typical properties and long-term evolution 
of core collapsed clusters is needed to make PE detection rates 
less uncertain. 

7.4. Implications 

Opening the gravitational-wave window to observe 
parabolic encounters of black holes in globular clusters of- 
fers a new possibility to constrain BH mass functions and GC 
models. Since PEs are very sensitive to the number of higher 
mass stellar BHs (Figs. 6 and 7), our results indicate that a reg- 
ular detection of PE events would provide accurate limits on 
the stellar BH mass distribution in GCs. Our analysis shows 
that this might be possible with AdLIGO if average GCs carry 
at least 500 BHs. 

7.4.1. Galactic Nuclei 

The analysis can be extended for other spherically symmet- 
ric systems using the scaling v oc N^n^^vJ^. Consider 
first galactic nuclei, hosting 2500 BHs and approximately all 
of these BHs have undergone mergers (Portegies Zwart & 
McMillan 2000). Galactic nuclei abundance in the universe is 
100 times less than for GCs (§ 4.2). Assuming that the virial 
velocity is a factor of \/T0 higher in galactic nuclei and that 
the CO mass function has the same distribution as in GCs, we 
get detection rates j/ gal = 5 2 x 1 /100 x 1 /VTO x i/ sc . However, 
the large number of BH mergers likely increases BH masses 
in galactic nuclei. For a uniform distribution (i.e. p = 0) 
of A^bh = 2500 between m m i n = 80 and m max = 100M Q we 
get 2500 2 /[500 x (100-80)/(100-5)] 2 x l/(100-s/l0) =1-8 
times the rates shown in the corresponding row of Tab. 1 
for GCs: 4. x lO^yr" 1 for InLIGO, 2.1 x 10" 3 for VIRGO, 
0.61 yr" 1 for AdLIGO, and 5.2 x lO^yr" 1 for LISA. These 
numbers should only be regarded as rough estimates, since 
they result from the direct application of simplified GC model 
assumptions to galactic nuclei. The calculation assumed un- 
corrected two-body interactions which does not hold for mo- 
tion in the potential of galactic centers (Rauch & Tremaine 
1996). 

7.4.2. Primordial Black Holes in Galaxies 

For a second example consider the GW detections from 
galactic haloes comprised of low-mass primordial BHs 
(PBHs) (see e.g. LIGO Scientific Collaboration 2005b and 
references therein). For a quick upper-limit estimate on the 
PE detection rate we repeat our analysis for GCs by changing 
the model parameters to describe galactic haloes. We assume 
NpBH = 10 11 PBHs within a maximum radius R = 5 kpc, a virial 
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velocity v v ; r = 220km/s, and a uniform distribution of masses 
between m m i n = O.25M0 and m max = O.95M0. The maximum 
distance of a matched filtering detection at a characteristic fre- 
quency of 0.9/Mmax = 4260 Hz with angular-averaged signal- 
to-noise ratio S/N = 5 is 0.33, 0.62, 4.5, 0, and 0.06 Mpc for 
InLIGO, VIRGO, AdLIGO, LISA, and NGLISA, respectively. 
The final result for the detection rate after the compilation of 
the full analysis described in the paper gives v = 1.5 x 10 -11 , 
1.6 x 10" 11 , 2.9 x 10" 10 , 2.6 x 10" 11 , and 1.4 x 10" 9 , respec- 
tively. These numbers are comparable to the total NS-NS PE 
rate in GCs. It is 9 orders of magnitudes smaller than the event 
rates estimates for PBH binary coalescence in one Milky Way 
sized galaxy (Ioka et al. 1998; LIGO Scientific Collaboration 
2005b). 

7.4.3. Unresolved Parabolic Encounter Background 

Another extension of the present analysis is to estimate the 
number of low S/N PE events, adding an unresolved astro- 
physical background to the GW detector noise budget similar 
to the unresolved WD background (e.g. Hils, Bender, & Web- 
bink 1990; Nelemans, Yungelson, & Portegies Zwart 2001; 
Benacquista, DeGoes, & Lunder 2004; Cornish & Crow- 
der 2005) and unresolved capture sources (Barack & Cutler 
2004b). Since PE rates are progressively larger for progres- 
sively smaller characteristic frequencies, fo, and since all PE 
waveforms extend to GW frequencies / < fo, PE background 
will be most substantial for space detectors, especially around 
the minimum frequency noise wall ( fmin = 10" 5 -10" 4 Hz). The 
total number of PE events within a distance D, can be obtained 
from (19) neglecting cosmology as 
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where / max is the maximum frequency for PEs avoiding a col- 
lision (32), which drops out to leading order if f m { a <C / max . 
ForD = lOGpc and / m ; n = 10" 5 Hz (lO^Hz) we get one event 
every \/R = 19 sec (88 sec), which corresponds to k = 5300 
(110) events for a 10~ 5 Hz (lO^Hz) frequency bin. If core 
contraction enhances PE rates by a factor of 14 4 (see § 7.3.3), 
we get k = 2 x 10 8 (4 x 10 6 ) events per frequency bin. 

These events will typically have a very low signal-to-noise 
ratio, e.g. for a single PE encounter with fo = 10" 5 Hz for 
mi=m 2 = 50M Q at d h = lOGpc we get (S/N)i - 8 x 10" 9 , 
and/o = lO^Hz yields (S/AOi - 2 x 10" 7 for LISA. Assuming 
that average unresolved PE noise accumulates proportional to 
\fk, we get a net amplitude of only (S/N) net = 10 -4 or 5 x 
10~ 4 for LISA even in the core contracted case for frequencies 
10~ 5 Hz or 10~ 4 Hz, respectively. Thus, we anticipate that the 
unresolved PE background adds a negligible amount to the 
LISA noise, and the unresolved PE background from stellar 
BH encounters will not be an issue for near-future extensions 
either. 
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APPROXIMATE ANALYTICAL FORMULAE 

The integrals given in Eq. (19) can be carried out analytically as a Taylor-expansion with respect to the small quantities 
^bh = m min/m m . dx and k N s = m N s/m max . The PE event rate is calculated in three parts 
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where ^f BH _ BH , ^fsH-NS' anc * ^hms-ns aie trie event rates of BH-BH, BH-NS, and NS-NS encounters. Substituting the mass 
dependence in Eq. (19) we get, 



v x (m u m 2 ) = 



v 2 (m\,m 2 ) = 



G 4 ' 3 N 2 C0 (m 1 + m 2 ) 4 / 3 



m 



3/2 



(47T) 2 /3 fl3 ( OT 7l+ OT -l)l/2 



,-2/3 
V100 ' 



2 1/3 G 2 / 3 JV ( 



CCvir 



(mi +m 2 ) 4/3 (m- l 1 +m 2 1 ) 1/2 ml> 2 f 



3/2 ,-4/3 



(A2) 



(A3) 



3tt 4 / 3 R 3 v " n ' "* z/ v " 1 ' '" 2 ' 100 
We approximate the gNs(fw) distribution with a Dirac-(5 function. Expanding the integrals in Taylor-series in k B h and k N s, we get 



^1, BH-BH 



= 2 -4/3. 
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2 5/6 



+ -j (1+2^+1.3644,,) 



y l ,NS-BH 
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1/3 
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(^bh) 
m C o 
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4/3 
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(A4) 



(w B h) 



1/3 



\ 2 

m NS \ l 



2 -5/2 ( «CO_A" 4/3 ( WnsV j 

\M Q J \M C J " 
where v\ is the Model-I event rate given by Eq. (17). 

Next we present ready-to-use formulas for calculating the detection rates of PEs. Assuming a constant density of GCs, no 
cosmological and no general relativistic corrections, the detection rate per logarithmic frequency bin becomes 

dv™ _ (2GM ) 19/3 2 „-3 v -iW3/Wo)\' A , rA7 . 



1, NS-NS " 



(A5) 
(A6) 
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where S/N is the minimum signal-to-noise ratio (which is set equal to 5 in our numerical results above), k= 85 3 / 2 2 1 / 3 7r 7//3 /6144 = 
2.323 is a constant coefficient, v v i r is the virial velocity (9), n gc is the average GC density in the universe, R gc is the typical radius 
of the GC, Nco is the number of COs in the GC, K(fo, A) and W(/o) are dimensionless terms, K(fo,X) depending on the CO mass 
distribution, gco, and W(fo) on the normalized GW energy spectrum F (f/fo) Eq. (5), and the detector spectral noise density, 
W>: 

K(J\).\)= / / d*id* 2 gco(*i)gco(*2) ( Xl+X2 y% » (A8) 




In terms of /o, K(fo,X) is constant for /o < /o,mox(2m m j n , A), decreases monotonically for larger fo and attains for /o > 
/o,max(2m max , A) (see Eq. [33] for the definition of /o,max(M,A), and m m ; n and m max are the minimum and maximum masses 
of the COs, respectively). In Eq. (A8), the integration variables x\ and X2 are the dimensionless masses of the COs, for which 
7Wmin/M Q < xi.2 < m max /M Q . For core collapsed clusters M has to be changed to m core , the typical mass of individual compo- 
nents in the core, R gc has to be changed to R COK , and v v i r to v core . These values should be set consistently with the core velocity 
dispersion and core radius which are input parameters for a globular cluster model. For given m components, the velocity disper- 
sion is then v m = (m/m C0K ) ' v core an d maximum radius from the cluster center is R m — (m/m cole ) 1//2 r? C ore- 

The total detection rate of parabolic encounters (again assuming a constant density of GCs and no cosmological and general 
relativistic corrections) is 



(2Gm V 9 / 3 /•/O.maxUMmax.-V 

k K ^e^ (5/A0 _3 n gc Ar 2 o/? -3 reV - ( l re / df f^W(f ) 3 K(f Q , A). (A10) 

C J/min/lO 



/•Amax(2M max ,A) 
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